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The  vector  components  of  several  baselines  shorter  than  100  meters  have 
been  measured  by  interferometric  processing  of  NAVSTAR  GPS  satellite  C/A-code 
broadcasts  to  an  accuracy  better  than  one  centimeter,  and  a  repeatability  on 
the  order  of  a  few  millimeters.  Interferometric  observables  were  extracted  by 
two  processing  schemes.  The  first  involved  broadband  cross  correlation  of  the 
modulated  GPS  .signals  received  by  the  antennas  that  define  the  unknown  base¬ 
lines.  The  second  approach  uses  GPS  receivers  at  each  antenna  to  reconstruct 
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-the  satellite  carrier  signal  from  which  phase  differences  can  be  measured.  In 
both  cases  the  interferometric  baseline  determinations  were  compared  in  a 
double-blind  test  with  a  conventional  survey  that  was  believed  to  be  accurate 
to  2  mm  (one-sigma)  and  3  arcseconds. 

Multipath  was  shown  to  be  a  significant  effect  only  to  the  extent  that  re¬ 
flecting  structures  might  block  the  visibility  of  useful  satellite  signals.  At 
the  present  time,  the  limit  on  extending  the  technique  to  precisely  measure 
baselines  exceeding  100  km  is  uncertainty  in  GPS  satellite  ephemerides.  For 
baselines  of  intermediate  length,  the  measurement  performance  is  limited  by  un¬ 
certainties  in  tropospheric  propagation  delays. 

The  results  of  the  experiments  plus  trends  toward  miniaturizing  GPS  equip¬ 
ment  complete  the  demonstration  that  short  baselines  can  be  measured  to  sub¬ 
centimeter  accuracy  using  highly  portable  equipment. 
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In  this  report  we  summarize  results  from  two  experiments  which 
were  designed  to  demonstrate  that 

(1)  Short  baselines  on  earth  can  be  measured  in  magnitude  and 
direction  within  one  (1)  centimeter  or  less  of  error,  using 
easily  portable  equipment,  and  that 

(2)  This  accuracy  can  be  achieved  by  interferometric  processing 
of  the  C/A-code  signals  currently  broadcast  by  NAVSTAR  GPS 
satellites . 

1.  EXPERIMENT  I:  Wideband  Signal  Processing 

Experiment  I  was  designed  to  demonstrate  the  extraction  of  inter¬ 
ferometric  observables  by  wideband  cross-correlation  of  the  GPS  signals 
received  at  each  end  of  a  baseline.  The  experiment  was  developed  when  it 
became  apparent  that  availability  of  the  GPS  receivers  required  for  Ex¬ 
periment  II  would  be  delayed  far  beyond  the  original  schedule. 

Experiment  I  was  planned  to  give  the  earliest  possible  credible 
demonstration  of  baseline  determination  using  GPS  signals.  This  was 
achieved  by  exploiting  resources  already  in  place  at  the  Westford,  Massa¬ 
chusetts  site  of  the  New  England  Radio  Observatory  Corporation  (NEROC) . 
These  facilities  included  the  MARK  III  VLBI  processor  which  was  used  to 
produce  "real-time"  data  representing  eight  lagged  correlation  products 
of  the  GPS  receiving  antenna  outputs  after  filtering  to  a  2  MHz  bandwidth 
encompassing  the  GPS  C/A  signal  at  1,574  Ghz. 


Other  objectives  of  the  experiment  were  to: 

(1)  Demonstrate  the  effectiveness  of  the  "MITES"  antenna  design 
(crossed-dipole  mounted  above  a  ground  plane)  which  pro¬ 
duces  a  nearly  omnidirectional  pattern  from  Zenith  to  with¬ 
in  10°  of  the  horizon. 

(2)  Demonstrate  the  effectiveness  of  a  novel  data-reduction  al¬ 
gorithm  which  approximates  a  "maximum  likelihood"  estimator 
of  the  baseline  vector  components. 

(3)  Explore  the  sensitivity  of  baseline  determinations  to 
oscillator  instability  and  multipath  to  the  extent  that 
these  corrupt  the  raw  data. 

The  Steinbrecher  Corporation  fabricated  the  "MITES"  antennas,  RF 
subsystems,  and  frequency  converters  to  condition  the  received  signals 
for  processing  by  the  MARK  III  unit.  Technical  direction  of  Experiment  I 
was  subcontracted  to  MIT  with  Professor  Charles  C.  Counselman  III  being 
the  Principal  Investigator.  Program  management  was  retained  by  Draper 
Laboratory.  Hie  final  report  prepared  by  Professor  Counselman  is  pro¬ 
vided  in  Enclosure  1 .  Unpublished  project  memoranda  are  provided  in  En¬ 
closure  2.  In  these  memoranda  Professor  Counselman  concludes  that,  "at 
least  when  the  satellite  geometry  is  good,  we  can  determine  all  three 
components  of  a  baseline  vector  that's  about  100  meters  long  within  about 
5  millimeters,  one  sigma,  by  MITES/GPS." 

2.  EXPERIMENT  II:  Narrowband  Signal  Processing 

The  concept  behind  Experiment  II  is  to  use  a  GPS  receiver  to  make 
precise  (albeit  ambiguous)  measurements  of  the  instantaneous  phase  of 
each  GPS  carrier  signal  in  the  narrow  bandwidth  of  the  carrier  tracking 
loop.  Carrier  phases  from  each  end  of  a  baseline  being  surveyed  were 
then  differenced  for  each  satellite  in  view  to  form  the  interferometric 
observables.  In  a  typical  experiment  the  observables  from  four  satel¬ 
lites  were  collected  every  1.92  s  during  a  period  of  about  one  hour. 

These  samples  were  recorded  on  a  computer  disk  and  then  processed  on  a 
minicomputer  ( TI  990/12)  to  estimate  the  baseline.  Approximately  30 
minutes  of  data  were  required  to  consistently  resolve  the  phase 
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ambiguities  that  affected  the  baseline  determination.  The  computer  pro¬ 
cessing  for  a  complete  baseline  determination  required  approximately  one 
hour  if  every  data  point  was  used;  experience  showed  that  good  results 
were  obtained  within  15-20  minutes  of  processing  time  if  only  every  10th 
sample  was  processed  (19.2  s  spacing). 

The  general  nature  of  the  results  is  that  baseline  measurements 
were  repeatable  at  the  level  of  1  mm  (standard  deviation)  in  the  horizon¬ 
tal  components  and  2  mm  in  the  vertical  component.  The  interferometric 
estimates  agreed  with  a  conventional  steel  tape  and  theodolite  survey  to 
within  3  millimeters  on  8  of  9  components  of  three  sample  baselines  and 
to  within  6  millimeters  on  the  ninth  component. 

Details  of  Experiment  II  and  the  analysis  of  its  results  are  re¬ 
ported  in  the  Symposium  paper  provided  as  Enclosure  3.  The  objectives  of 
Experiment  II,  which  were  achieved,  were: 

(1)  To  demonstrate  the  feasibility  of  sub-centimeter  surveying 
of  short  baselines  using  easily  portable,  and  relatively 
affordable  observing  equipment. 

(2)  To  assess  the  limits  of  baseline  measurement  accuracy  and 
repeatability,  and  to  identify  the  dominant  error  sources. 
With  respect  to  the  first  objective  we  demonstrated: 

(a)  The  independent  measurement  of  three  legs  of  a  tri¬ 
angle  with  vector  misclosure  of  (-2.2,  -2.2,  -0.3) 
mm  in  the  coordinate  directions  (North,  East,  Up). 

(b)  Repeatability  of  one  baseline  measurement  over  four 
independent  trials  at  the  1-sigma  level  of  (0.76, 
1.60,  2.40)  mm  in  the  (N,  E,  U)  coordinates. 

(c)  The  largest  disagreement  between  any  interferometri- 
cally  measured  baseline  component  and  the  corre¬ 
sponding  component  measured  by  a  conventional  survey 
was  5.7  mm.  Five  of  the  nine  baseline  components 
were  measured  with  less  than  a  0.4  mm  discrepancy. 
The  nominal  uncertainty  of  the  conventional  survey 
was  given  as  2  mm  (one-sigma). 

Efforts  are  currently  underway  under  sponsorship  of  the  Defense 
Department  Advanced  Research  Project  Agency  to  miniaturize  GPS  receivers 


to  approximately  the  size  of  a  ciqarette  pack..  It  is  easy  to  visualize 
that  by  the  late  1980s  short-baseline  survey  equipment  will  consist  of  a 
small  portable  antenna,  an  integrated  RF/GPS  receiver  assembly,  a  crystal 
oscillator  (with  short-term  stability  exceeding  one  part  in  lO1^  over 
100  sec),  a  digitizer/recorder  or  digitizer/modem  unit,  and  power  sup¬ 
plies.  The  electronics  would  be  packaged  in  a  small  suitcase.  Carrier 
phase  measurements  from  the  four  satellites  at  a  rate  of  3  samples  per 
satellite  per  minute  over  one  hour  of  observations  should  be  adequate  to 
achieve  sub-centimeter  accuracy.  An  inexpensive  storage  medium  or  commu¬ 
nication  link  capable  of  handlinq  a  total  of  about  1000  16-bit  words  for 
each  baseline  is  all  that  is  required  to  service  data  collection  require¬ 
ments  . 

With  respect  to  the  second  objective,  we  observed  that  the  post- 
fit  residual  error  for  measurement  of  a  null  (zero-length)  baseline  was 
0.96  mm,  compared  to  0.55  mm,  which  was  predicted  from  the  tracking  loop 
bandwidth  and  the  estimated  signal-to-source  ratio.  This  is  satisfactory 
agreement  given  the  uncertainty  in  the  signal-to-noise  ratio  measurement. 

The  post-fit  residual  error  for  baseline  measurements  using  carri¬ 
er  phase  observations  was  about  2  mm.  Thus  the  residual  was  approximate¬ 
ly  twice  that  of  the  noise  level  of  the  observing  system.  The  increase 
is  associated  to  a  great  extent  with  periodic  variations  which  affect  all 
satellite  observations  in  the  same  fashion.  The  period  of  these  varia¬ 
tions  is  about  forty  minutes  (see  Enclosure  3,  Figure  7).  The  short-term 
scatter  around  the  periodic  variation  is  comparable  to  the  noise-only  re¬ 
siduals  shown  in  Figure  5  of  Enclosure  3).  We  attribute  the  periodic 
components  to  imperfect  modeling  of  the  differential  expansion  of  the 
cables  connecting  each  antenna  to  the  data  collection  system,  but  this 
has  not  been  conclusively  established.  If  true,  it  would  suggest  that 
field  survey  equipment  should  use  the  smallest  possible  length  of  cable 
to  connect  the  GPS  antennas  to  the  GPS  receivers.  At  any  rate,  the  2  mm 
residual  is  well  within  our  projected  expectation. 

Other  factors  that  once  were  thought  to  limit  the  performance  in¬ 
clude  oscillator  drift  and  multipath.  However,  the  combination  of  high 
quality  quartz  crystal  oscillators  and  the  use  of  a  sequential  GPS  re¬ 
ceiver  which  can  sample  the  carrier  phase  to  each  satellite  in  view 


faster  than  once  per  second  combine  to  make  oscillator  drift  a  negligible 
problem.  For  example,  if  the  Allen  variance  of  the  oscillator  over  1 
second  is 


then  the  contribution  of  clock  drift  over  1  second  to  the  differential 
measurement  of  carrier  phase  is  /2  x  10-12f  cycles  rms,  which  is  less 
than  0.003  cycles  (0.6  mm)  at  the  GPS  LI  frequency.  At  this  level,  clock 
drift  would  be  about  five  times  less  significant  than  thermal  noise.  If 
one  wished  to  economize  by  using  a  lower  quality  oscillator,  then  some  of 
the  increased  frequency  jitter  could  be  offset  by  faster  cycling  among 
GPS  satellites.  For  example,  one  commercially  available  GPS  receiver  cy¬ 
cles  among  four  satellites  every  20  ms. 

The  primary  effect  of  multipath  will  probably  be  blockage  of  sat¬ 
ellite  visibility,  rather  than  distortion  of  signal  phase  measurements. 
Enclosure  4  presents  some  calculations  based  on  an  analysis  by  Counselman 
which  predict  the  carrier  phase  fluctuations  attributable  to  a  reflecting 
surface  mounted  in  the  vicinity  of  a  "MITES"  antenna.  The  specific  con¬ 
figuration  considered  in  that  analysis  was  set  up  in  an  experiment  to 
check  the  theoretical  predictions.  The  results  (compare  Figure  10  of  En¬ 
closure  4  with  Enclosure  3)  were  qualitatively  consistent  for  all  satel¬ 
lite  observations  and  within  10%  error  for  the  satellite  that  suffered 
the  greatest  multipath  effects.  We  conclude  that  the  effects  of  multi- 
path  will  average  out  of  any  series  of  GPS  carrier  phase  measurements 
taken  over  a  half  hour  or  longer  interval;  moreover,  any  reflecting  sur¬ 
faces  that  are  greater  than  1000  feet  from  the  antenna  (one  C/A-code 
"chip")  will  produce  multipath  that  is  reduced  below  thermal  noise  level 
by  the  GPS  signal  processing.*  Therefore  multipath  should  not  be  a  sig¬ 
nificant  problem  for  well  engineered  observing  systems  with  good  siting 
of  the  receiving  antennas.  Exceptions  might  occur  if  one  had  to  measure 
a  baseline  extending  from  a  large  structure,  or  in  the  vicinity  of  metal 
towers. 

*One  could  achieve  immunity  to  closer  sources  of  multipath  by  processing 
GPS  P-code  signals,  for  which  the  "chip"  width  is  equivalent  to  a  100 
foot  span.  The  costs  of  a  P-code  receiver  would  be  higher  than  for  a 
C/A-code  receiver. 
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It  appears  that  the  fundamental  limits  of  baseline  measurement  ac¬ 
curacy  will  be  thermal  noise  in  the  GPS  receiver,  ephemeris  errors  in 
location  of  GPS  satellites  and  excess  propagation  delay.  If  the  uncer¬ 
tainty  in  satellite  location  is  Ap  (meters),  then  the  resuting  uncertain¬ 
ty  A£  in  measuring  a  baseline  whose  length  is  £  is  approximately 

A£  =  -L4L- 
26  x  106 

Baseline  measurement  is  most  sensitive  to  "along  track"  errors  in  satel¬ 
lite  locations  and  less  sensitive  to  radial  errors.  For  current  GPS 
broadcast  ephemerides  that  are  no  more  than  1  day  old,  the  "along  track" 
error  (one-sigma)  is  about  20  meters;  hence,  A £/£  a  10~®.  Clearly,  im¬ 
provement  in  ephemeris  modeling  will  be  needed  to  preserve  centimeter 
level  accuracy  at  distances  in  excess  of  about  10  km. 

Propagation  delays  sufferred  by  GPS  signals  as  they  traverse  the 
ionosphere  and  troposphere  en  route  from  satellites  to  ground  terminals 
are  highly  correlated  at  antennas  separated  by  short  distances.  Thus, 
the  effects  of  these  delays  are  essentially  cancelled  by  phase-difference 
processing.  As  baseline  lengths  increase,  the  delays  tend  to  decorre¬ 
late.  The  impact  of  ionospheric  errors  may  be  made  negligible  by  dual¬ 
frequency  interferometry  on  GPS  Lj  and  L2  frequencies. 

Because  uncertainties  in  modeling  tropospheric  propagation  delays 
vary  from  about  10  cm  rms  at  90°  elevation  to  about  50  cm  rms  at  20°  ele¬ 
vation,*  and  because  there  can  occasionally  be  larqe  variations  in  pres¬ 
sure,  temperature,  and  humidity  over  short  distances,  tropospheric  errors 
could  introduce  1-cm  baseline-measurement  errors  with  baselines  as  short 
as  1  km.  However,  when  weather  conditions  are  quite  uniform  and  dry  at 
both  ends  of  the  baseline  (and  along  the  propagation  paths),  and  when 
meteorological  parameters  are  measured  and  incorporated  in  the  propaga¬ 
tion  model,  baseline-measurement  errors  presumably  can  be  maintained  to 
less  than  1  cm  f or  baselines  much  longer  than  1  km.  Under  good  weather 
conditions,  the  baseline-measurement  errrors  due  to  tropospheric  effects 

•Using  measured  values  of  the  surface  index  of  refraction. 


can  be  maintained  within  about  10  cm  for  arbitrarily  long  baselines.** 
For  long  baselines,  the  1  part-per  million  errors  due  to  GPS  ephemeris 
uncertainties  become  dominant. 


3.  CONCLUSIONS: 

The  two  experiments  carried  out  under  this  program  have  demon¬ 
strated  the  use  of  GPS  signals  for  the  measurement  of  short  baselines 
with  accuracies  of  better  than  1  cm  rms.  With  the  continued  evalution  of 
electronic  technology,  the  cost  of  interferometry  eguipment  for  such 
measurements  is  expected  to  drop  dramatically  in  the  near  future,  so  that 
both  of  the  described  techniques  will  be  seen  to  be  economically  attract¬ 
ive  as  compared  with  alternative  surveying  techniques. 

The  second  technique,  utilizing  GPS  C/A-code  receivers,  offers  the 
advantage  of  providing  absolute,  as  well  as  relative,  positioning  data. 
Moreover,  with  the  projected  application  of  such  receivers  to  automobiles 
and  other  mass  markets,  the  associated  equipment  costs  are  expected  to 
become  particularly  attractive. 


* *As  indicated  by  typical  VLBI  results  reported  since  1978.  For  example, 
see  "Analysis  of  Lunar  Laser  Ranging  Data  and  Performance  and  Analysis  of 
VLBI  Observations  for  Geodetic  Purposes,"  AFGL-TR-81-004,  C.C.  Counsel- 
man  III,  et  al. ,  27  November  1980. 
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Under  this  subcontract  m.I.T.  planned,  performed,  and 
analyzed  radio  interferometric  observations  of  N’AVSTAR/GPS 
satellites. 

In  a  25-page  memorandum  dated  June  11,  1980,  to  R.  L. 
Greenspan  from  C.  C.  Counselman  III  and  S.  A.  Gourevitch,  we 
described  the  observational  data  that  we  expected  to  acarire,  and 
some  alaorithms  that  we  intended  to  use  for  analysis  of  these 
data . 

In  a  second  15-page  memorandum  dated  August  13,  1980,  to 
R.  L.  Greenspan  from  C.  C.  Counselman  III  and  S.  A.  Gourevitch, 
we  further  described  the  data-analysis  algorithms. 

In  a  19-paae  memorandum  dated  September  30,  1980,  to  R.  L. 
Greenspan  from  C.  C.  Counselman  III  and  S.  A.  Gourevitch,  we 
reported  on  a  computer  simulation  of  interferometric  observations 
in  which  multipath  effects  were  included. 

Actual  observations  of  the  GPS  satellites  were  performed  at 
Haystack  Observatory  in  December,  1980,  and  subsequently  were 
analyzed  to  determine  interferometer  baseline  vectors.  These 
experiments  and  the  initial  analysis  results  were  reported  at  the 
1981  Spring  meeting  of  the  American  Geophysical  Union  by  C.  C. 
Counselman  III.  An  abstract  of  this  presentation,  entitled 
"Accuracy  of  baseline  determinations  by  MITES  assessed  by 
comparison  with  tape,  theodolite,  and  geodimeter  measurements," 
was  published  in  Eos  (Trans.  AGU),  vol.  62,  p.  260,  on  April  28, 
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Another  presentation  was  made  by  C.  C.  Counselman  III  to  the 
1981  IEEE  International  Geoscience  and  Remote  Sensina  Symposium, 
ana  appears  on  pp.  219-224  of  vol.  1  of  the  Symposium  Digest . 

Our  data-analysis  algorithm,  multipath  simulation  results, 
and  a  few  actual  experimental  results  were  also  published  in  a 
journal  article  by  C.  C.  Counselman  III  and  S.  A.  Gourevitcn, 
entitled  "Miniature  interferometer  terminals  for  earth  surveying: 
ambiguity  and  multipath  with  global  positioning  system," 
appearing  in  IEEE  Transactions  on  Geoscience  and  Remote  Sensing, 
vol.  GE-19,  pp.  244-252,  October  1981. 

The  results  of  a  comparison  of  our  GPS  radio  interferometry 
determinations  of  a  triangle  of  baselines  with  an  independent 
determination  by  conventional  survey  methods  were  given  in  detail 
in  a  9-page  memorandum  to  R.  L.  Greenspan  from  C.  C.  Counselman 
III  dated  September  30,  1981. 

Copies  of  all  of  the  documents  mentioned  have  already  been 
provided  to  Dr.  R.  L.  Greenspan,  who  was  the  technical  monitor  of 
this  subcontract. 

No  inventions  were  made  in  work  under  this  subcontract. 
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DEPARTMENT  OF  EARTH  AND  PLANETARY  SCIENCES 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 

CAMBRIDGE.  MASSACHUSETTS  021S9 

54-626 

June  11,  1980 

Memorandum 

To:  R.  L.  Greenspan 

From:  C.  C.  Counselman  and  S.  A.  Gourevitch 

Subject:  Data  Processing  Algorithms  for  Experiment  I — 

With  the  Haystack  Mark  III  Back  End  and  Correlator 

1.  Introduction 

In  this  memorandum  we  describe  the  data  that  we  expect 
to  acquire,  and  algorithms  that  we  expect  to  use  for  the 
analysis  of  this  data,  in  the  "Experiment  I"  described  in 
Section  4.1  of  your  DL  Intralab  Memo  no.  94100-01A.  In  this 
experiment,  signals  received  from  GPS  satellites  by  antennas 
in  the  Haystack  Observatory  parking  lot  will  be  processed 
by  Haystack's  Mark  III  VLBI  correlator. 

2.  Block  Diagram 

Figure  1  shows  a  block  diagram  of  some  of  the  equipment 
to  be  used  in  the  experiment.  The  GPS  Lj  signals  received  by 
one  antenna,  at  nominal  signal  carrier  frequency  1575.42  MHz, 
with  C/A  and  P  code  modulation,  and  with  Doppler  shifts,  will 
be  converted  in  two  stages  of  mixing  to  a  "video"  frequency 
band,  nominally  0-2  MHz,  and  sampled  once  per  0.25  usee  to 
yield  a  series  of  samples  x^  at  uniformly-space  times  t^. 

From  another  antenna  in  a  similar  way  we  obtain  a  series 
of  samples  y^  at  essentially  the  same  times.  Before  being 
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sampled,  the  signals  are  clipped  so  that  only  their  signs 
are  preserved?  that  is,  the  sample  values  are  ±1  and  can 
be  represented  by  just  one  bit  per  sample. 

A  single  video  passband  only  2  MHz  wide  can  contain 
most  of  the  C/A  modulation  power  but  only  about  20%  of  the 
P  power.  In  order  to  capture  more — in  fact  most— of  the 
P  power,  we  will  use  several  2-MHz-wide  video  channels  in 
parallel,  each  with  a  different  second  local  oscillator 
frequency.  These  frequencies  will  be  spaced  by  integer 
multiples  of  2  MHz.  To  keep  the  present  discussion  simple, 
however,  we  will  pretend  that  there  is  only  one  channel,  and 
we  will  ignore  the  problem  of  how  to  combine  the  data  from 
different  channels.  (This  is  an  old  problem,  with  a  known 
and  straightforward  solution.) 

Not  shown  in  Figure  1  is  the  system  with  which  the 
sampled,  time- tagged,  one-bit  data  are  formatted,  recorded  on 
magnetic  tape,  and  later  reproduced  for  processing  by  the 
Mark  III  correlator. 

3.  Mark  III  Correlator 

The  Mark  III  correlator  computes  eight  values  of  the 
crosscorrelation  function  between  the  two  sampled  signals, 
for  eight  discrete  values  of  delay  offset  between  the  sample 
streams.  It  computes,  once  every  2  seconds  of  time  t,  the 
2-second  averages 

r£(^J_^2Kiyi+Jl  ;  1  =  ~4'  ”3'  ~2'  +3 

^successive) 

where  N  =  8x10 5  is  the  number  of  products  summed;  the 
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time  span  of  2  seconds  {=  8x10s  x  0.25  usee)  is  centered  on 
the  time  t.  For  simplicity,  we  are  deliberately  suppressing 
the  book-keeping  details,  as  well  as  some  calculus  that  shows 
that,  with  sufficient  accuracy,  r^(t)  computed  from  one-bit 
samples  is  almost  equivalent  to  what  you  would  get  had  you 
had  an  infinite  number  of  bits.  (The  only  significant  lack 
of  equivalence  is  that,  with  one  bit,  the  signal-to-noise  ratio 
is  reduced  by  a  factor  of  2/ir.) 


4.  Expected  Form  of  r 


Ignoring  some  of  the  noise"  effects,  the  expected  func¬ 
tional  form  of  r^  is 
K 

r£  "  2  *k  *  MTk“£6)  *  cosI“cTk  “*LO  ~(“c"WLOH65‘' 


(Eqn.  1) 


where  k  is  the  index  number  of  a  satellite  (several  may  be 
visible  simultaneously) ;  A^  is  the  ratio  of  the  peak  (at  the 
carrier  frequency)  received  power  spectral  density  of  the  kth 
satellite  signal  to  the  background  noise  power  spectral  density; 
A (x)  is  a  (triangle)  function  equal  to  the  autocorrelation 
function  of  the  C/A  modulation  (we  ignore  for  now  the  P  modu¬ 
lation)  ; 


A  (x) 


|l  -1.023  MHz  •  |x|  | 


whichever  is  greater; 


(Eqn.  2) 


is  the  group  delay  difference  between  the  kth  satellite 
signals  at  the  x  and  y  samplers;  6  is  the  inter-sample  time 
interval,  equal  to  0.25  us;  u  is  the  transmitting  carrier 

C 


6 


frequency,  in  radians  per  second,  assumed  the  same  for  all 

transmitters;  4>LQ  is  the  difference  between  the  phases  of  the 

two  different  receivers'  local  oscillators;  <tsLQ  is  the  (radian 

frequency  of  the  local  oscillators,  assumed  the  same  for  both 

receivers;  the  difference  u  -<oT_  is  equal  to  2ir  x  10 6  radians 

c  LO 

second.  The  time-variations  of  the  various  terms  in  Eqn.  1 
have  not  been  shown  explicitly,  but  need  to  be  kept  in  mind. 
We  discuss  each  term  briefly. 

The  "amplitude"  A^  will  vary  slowly  due  to  the  variation 
of  the  receiving  (and  to  a  lesser  extent,  the  transmitting) 
antenna  gain  with  elevation  angle,  coupled  with  the  motion 
of  the  satellite  in  the  sky.  The  variation  will  be  by  less 
than  a  factor  of  about  4  for  elevation  angles  from  10°  to 
90°.  For  elevation  angles  above  about  20°,  our  antenna  has 
at  least  unity  gain  and  the  received  peak  C/A  power  spectral 
density  should  be  at  least  10*"22  watt/Hz.  The  noise  power 
spectral  density  with  a  5  db  noise  figure  will  be  about 
1.25xl0“20  watt/Hz.  Therefore,  we  can  expect  A^  -  8xl0~s 
for  the  C/A  modulation.  For  the  P  modulation,  which  we  are 
temporarily  ignoring,  A^  =  4X10-11. 

The  group  delay  combines  the  free-space  propagation 
delay  and  delays  within  the  receivers  and  cables.  All  ingred: 
ents  of  t,  are  differenced  between  the  x  and  y  systems.  The 
greatest  contributor  to  and  to  its  time-derivative  is 
likely  to  be  the  free-space  delay.  The  magnitude  of  this 
delay  contribution  will  be  of  the  order  of  the  baseline  lengt) 
measured  in  light-time;  that  is,  for  example,  about  300  ns  for 
a  100-meter  baseline.  The  delay  rate,  dxt/dt,  will,  be  of  the 


5 


order  of  the  baseline  length  times  the  apparent  satellite 
orbital  angular  velocity,  about  l^xlO-4  radian/second  maxi¬ 
mum.  The  peak  delay  rate  for  a  100-meter  baseline  is  thus 
about  4x10  11  sec/sec,  and  the  corresponding  "fringe  rate", 
<ocTjc,  about  0.4  radian/second.  The  corresponding  minimum 
fringe  rotation  period  is  about  15  seconds,  and  the  maximum 
amount' of  rotation  during  one  2-second  integration  is  7O.8 
radian  *  47  degrees,  which  is  just  tolerable:  the  coherence 
loss  due  to  the  rotation  during  the  integration  is  about  3%. 

The  difference  between  the  x  and  y  receiver  local-oscil¬ 
lator  phases  (and  any  other  instrumental  contributions  to 
the  phase  delay  in  excess  of  the  group  delay  t^) ,  represented 
by  <f>LQ,  is  expected  to  be  a  slowly  varying  function  tof  time, 
changing  by  less  than  does  over  any  relevant  time  interval.  All 

of  the  instrumental  phase  drift  is,  by  definition,  lumped 
into  the  time-variation  of  4^.  Thus,  is  a  constant  in 
Eqn.  1. 


It  is  also  permissible  to  regard  u  as  a  constant.  The 
actual  satellite  carrier  frequencies  are  unlikely  to  deviate 
from  the  nominal  value  of  1,575.42  MHz  by  more  than  1  kHz. 

The  maximum  magnitudes  of  and  of  ii  will  be  t  1  ysec. 
Therefore  the  argument  of  the  cosine  function  in  Eqn.  1  cannot 
be  in  error  by  more  than  about  10“ 5  cycles  of  phase  by  virtue  of 
the  error  in  the  assumption  of  constant  coc  equal  to  the  nomi¬ 
nal  value. 

» 

Since  the  magnitude  of  will  be  limited  by  the  base¬ 
line  length,  for  example  to  about  500  ns  for  a  500-foot  base¬ 
line,  the  peak  of  the  triangle  function  will  occur  in 

this  example  for  ,£|<  2.  (.Recall  that  6  =  250  ns.)  The 
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half-width  of  the  base  of  the  triangle  is  approximately 
1000  ns;  therefore,  even  with  an  extreme  value  of  most 
of  the  area  of  the  triangle,  corresponding  to  most  of  the 
"weight"  of  information  in  the  r^ ,  will  be  contained  within 
the  range' -4<£<+3  for  which  r^  is  computed  by  the  Mark  III 
correlator.  If  we  should  wish  to  use  a  baseline  much 
longer  than  500  feet,  we  would  have  two  options  for  extend¬ 
ing  the  range  of  S.6.  The  first  option  is  to  increase  6. 

This  can  be  done  in  factors  of  2  (to  a  maximum  5*4  ps)  , 
with  a  concomitant  reduction  of  the  bandwidth  of  a  video 
channel.  The  second,  more  attractive,  option  is  to  make 
repeated  passes  of  the  tape-recorded  data  through  the 
correlator.  On  each  pass,  the  r^  are  computed  for  eight 
values  of  £,  but  a  different  set  of  eight  values  can  be 
obtained  on  different  passes.  (However,  on  each  pass  the 
eight  values  must  be  successive  integer  values.) 

The  form  of  r^  as  a  function  of  £,  for  a  single  value 
of  k,  for  equal  to  zero,  and  for  two  different  values  of 
4  _,  equal  to  0  and  r/2,  is  shown  in  Figure  2.  With  the 
summation  over  k  and  with  a  variety  of  values  of  t^, 
obviously  the  picture  can  become  very  messy.  With  a  base¬ 
line  length  under  500  feet,  all  the  triangle  functions  overl 
and  with  the  various  t^'s  differing  by  "random"  amounts,  the 
relative  phases  of  the  cosines  are  random.  Extracting  order 
from  this  overlapping  mess  is  the  job  of  our  data  analysis 
algorithm. 


5.  Noise 


The  extraction  must  be  done  in  the  presence  of  signifi¬ 
cant  amounts  of  background  noise,  which  really  i£  random. 
Again  considering  that  only  one  satellite  is  present,  to 
simplify  the  discussion,  we  can  show  that  the  signal-to- 
noise  ratio  at  the  expected  peak  of  r^-  (that  is,  with  £6  = 
and  with  Wj «=  -4>LQ)  is 

SNR  -  |  •  N*^  •  Ag  (Egn.  3) 

where  the  factor  (2/ir)  represents  the  one-bit  clipping 
loss,  N  is  the  number  of  correlated  pairs  of  samples  (equal 
to  8x10*  for  a  single  2-MHz  video  channel  and  a  2-second 
integration  time) ,  and  Ag  is  the  ratio  of  the  total ^signal 
power  contained  in  the  channel  bandwidth  to  the  noise  power 
in  this  band.  For  the  C/A  modulation  and  a  2-MHz  bandwidth, 
Ag  —  0.5  times  the  peak  signal-to-noise  power  spectral 
density  ratio  represented  by  A^  in  Eqn.  (1).  Thus,  the 
expected  single-satellite  C/A  SNR  will  be  -7  for  the  band¬ 
width  and  integration  time  given  above.  (For  the  P  signal, 
the  combination  of  the  20- times- lower  peak  power  spectral 
density  and  the  10-times-higher  bandwidth  reduces  the  SNR 
by  a.  multiplicative  factor  of  20**^  •  10*5  *  O'.  16,  if  the 
number  of  video  channels  used  is  increased  from  1  to  10 
in  order  to  capture  the  same  fraction  of  the  modulation 

power.  It  should  be  noted,  however,  that  the  best  SNR  is 

* 

obtained  with  less  than  20  MHz  bandwidth  for  the  P  modulation 
and  less  than  2  MHz  for  the  C/A.) 

The  standard  deviation  of  the  error  in  the  estimate  of 
the  "fringe  phase"  (to  t.  )  ,  in  radians,  is  just  the 
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reciprocal  of  the  SNR,  in  the  single-satellite,  noise- 
limited  case. 

The  standard  deviation  of  the  error  in  the  estimate 
of  the  group  delay,  t^,  is  approximately  independent  of  the 
video  channel  bandwidth  or  the  number  of  channels  combined, 
as  long  as  the  combined  bandwidth  exceeds  about  1.6f^r 
where  ffa  =  1.023  MHz  for  the  C/A,  and  ffe  =  10.23  MHz  for 
the  P  modulation. .  We  find  that 


"t  =  T  <  (Eqn.  4) 

where  T  is  the  integration  time  used.  Thus,  for  the  C/A 
signal  alone. 


VC/A)  =  100  ns  •  T^conds  , 


and  for  the  P  signal  alone, 

nT(P|  =  62  ns  •  ^seconds  ' 


(EqnV  5) 


(Eqn.  6) 


where  T-seconas  ^-s  the  integration  time  in  seconds.  Accord¬ 
ing  to  these  formulae,  there  is  not  a  very  great  net 
advantage  to  using  the  wider-bandwidth  P  code  modulation, 
as  opposed  to  the  C/A.  The  advantage  of  the  10  times 
wider  bandwidth  is  largely  canceled  by  the  disadvantage  of 
the  lower  signal-to-noise  ratio.  In  practice,  one  expects, 
the  wider-bandwidth  approach  would  pay  off  because  one 
would  not  be  trying  to  "split  a  chip"  so  finely.  Therefore, 
one  would  not  be  as  susceptible  to  systematic  errors  such 
as  those  due  to  multipath. 


6.  Separating  Satellites 

From  the  preceding  discussion  it  should  be  apparent 
that  the  worse  problem  in  the  data  analysis  for  the  planned 
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experiment  is  likely  not  to  be  background  noise,  but 
rather  will  be  interference  between  satellites.  The 
interference  problem  forces  us  to  use  a  rather  complicated 
algorithm  in  order  to  extract  the  best  estimate  of  -the 
baseline  vector,  etc.,  from  the  r^‘  data.  We  are  consid¬ 
ering  two  different,  basic,  approaches which  we  now  outline 

7.  Traditional  Approach 

The  traditional  approach  to  the  analysis  of  VLBI 
correlation  data  (as  described,  for  example,  by  A.  R. 

Whitney  et  al. ,  A  very-long-baseline  interferometer  system 
for  geodetic  applications ,  Radio  Science,  11,  421-432,  1976) 
makes  use  of  a  priori  information  on  the  baseline  vector, 
the  positions  and  motions  of  the  sources,  etc.,  in  order  to 
"counter-rotate"  the  "fringes"  and  to  shift  the  y  signal  in 
group  delay  relative  to  the  x  signal.  If  we  have  accurate 
information  on  the  baseline  and  the  satellite  ephemerides 
in  our  experiment,  then  in  our  data  analysis  we  can  make 
the  fringes  "stand  still"  for  any  one  satellite  at  a  time. 
That  is,  the  contribution  of  one  satellite  to  the  sum  over 
k  in  Eqn.  1  can  be,  in  effect,  held  nearly  constant  while 
the  contributions  from  the  other  satellites  "rotate". 
(Actually,  since  Eqn.  1  contains  only  the  real-valued  cosine 
function  and  not  the  complex  exponential  usually  seen, 
one  might  prefer  to  say  that  our  fringes  don’t  rotate, 
they  just  travel  along  the  l  axis  in  Figure  2,  as  the 
fringe  phase  advances.  We  shall  use  the  traditional  jargon.) 
If  the.  fringe  rate  difference  between  the  satellite  chosen  to 
stand  still  and  another,  potentially  interfering,  satellite 


is  sufficiently  high,  the  interfering  fringes  will  wash 
out,  or  average  to  zero,  well  enough  over. the  integration 
time. 

This  approach  has  two  drawbacks  in  the  context  of 
our  experiment.  First,  it  requires  fairly  accurate  infor¬ 
mation  on  the  baseline  vector;  this  information  could  be 
expensive  to'  get,  especially  if  we  are  going  to  move  our 
antennas  many  times.  For  the  a  priori  information  to  be 
useful  in  discriminating  between  satellites,  its  fractional 
uncertainty  must  be  small.  The  effective  interference 
attenuation  factor  is  approximately  equal  to  the  fractional 
(dimensionless)  uncertainty  in  the  a  priori  prediction  of 
the  fringe  rate  difference.  Probably  we  would  needy  to 
have  the  magnitude  of  the  a  priori  uncertainty  of  our 
baseline  vector  be  less  than  a  few  percent  of  the  baseline 
length.  Also,  since  the  fringe  rate  difference  between 
two  satellites  would  occasionally  be  zero  or  nearly  zero, 
there  would  be  times  when  the  two  satellites'  fringes  were 
difficult  or  impossible  to  distinguish  within  reasonable 
integration  times.  (Some  advantage  could  be  taken  of 
group  delay  resolution,  since  the  group  delay  difference 
between  two  satellites  will  be  relatively  large  when  the 
delay  rate  difference  is  small.  However,  this  advantage 
will  be  marginal  when  the  group  delay  difference,  which  can 
never  exceed  twice  the  baseline  length  in  light-time,  is 
less  than  the  C/A  code  "chip"  width,  of  about  1,000  ns.) 
Finally,  the  requirement  of  substantial  a  priori  inputs 
might  be  considered  philosophically  awkward,  in  the  sense 
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of  "What  good  is  your  system  if  it  needs  to  know  the  answer 
in  advance?" 

The  second  drawback,  am  extension  of  the  first,  is 
that,  since  the  coherent  integration  time  required. to 
filter  out  the  interference  is  inversely  proportional  to  the 
fringe  rate  difference,  the  time  required  is  inversely 
proportional  to  the  baseline  length  and  becomes  unreasonably 
large  for  baselines  shorter  than  about  100  feet.  With 
100-foot  length,  the  maximum  possible  fringe  rate  is  about 
0.02  Hz,  and  a  substantial  fraction  of  the  time  there  will 
be  two  satellites  with  fringe  rates  differing  by  as  little 
as  one-tenth  the  maximum,  or  0.002  Hz.  The  beat  period  in 
this  case  is  500  seconds,  and  an  integration  period  several 
times  longer  than  this  beat  period  will  be  required  for 
effective  separation  of  the  two  satellites. 

We  conclude  that  we  should  take  an  approach 
different  from  the  traditional  one. 

8.  C4  Method 

We  suggest  an  approach  based  upon  direct  maximization 
of  the  cross-correlation  of  the  observed  cross-correlations 
r^(t)  with  a  theoretical  model.  We  dub  this  approach  the 
C*  method.  The  idea  is  elementary.  Using  Eqn.  (1)  ,  a  set 
of  theoretical  values  r^(t)  may  be  calculated  from  any  set 
of  assumed,  trial  values  of  the  unknown  parameters:  the 
three  (constant)  components  of  the  baseline  vector  in  an 
Earth-fixed  coordinate  system,  and  the  (also  assumed  constant) 
instrumental  group  delay.  These  four  unknowns  are  fundamental, 
and  enter  Eqn.  (1)  through  t^.  As  a  first  approximation. 
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the  amplitudes  may  be  treated  as  constants,  or  at  least 
as  knowns,  calculated  from  the  elevation  angles  and  the 
nominal  antenna  gain  patterns.  The  quantities  6,  coc,  and 
w^Q.are  constants  known  with  more  than  sufficient  accuracy. 

The  phase  4>L0  is  unknown  but,  because  it  is  not  a  function 
of  k,  it  may  be  suppressed  by  replacing  the  cosine  function 
in  r^  (t)  with  the  complex  exponential ,  and  then  maximizing  the 
magnitude  of  the  complex  correlation  between  r^(t)  and 
r.£(t).  It  is  necessary,  in  order  to  suppress  <j>LQ  exactly, 
that  be  constant  over  the  time  interval  of  the  "Ch" 

LiU 

operation.  To  summarize,  we  shall  compute 


T,  =  X.  (B,  X  ,  t)  , 


k  =  1**K, 


(Eqn.  7) 


where  x^  is  a  function  of  the  trial  value,  or  "guess",  B  of 

the  baseline  vector  S,  the  trial  value  x1,,  of  the  instrumental 

group  delay  x  ,  and  the  known  time,  t.  Then  we  compute  the 
o 


complex  correlation 


C  =  C(§,  J  ^  <*t)  (t)  cdt 


(Eqn.  8) 


where  the  time  span  of  the  integration  (really  a  sum  over 

discrete  2-second  averages)  can  be  quite  long,  perhaps  the 

duration  of  the  entire  satellite  "pass,"  and  where 
K 

r£(t)  =2  ^(t).  •  A[xk(t)  -16]  •  exp[ja>cxk  -j  (“c-wL0H$)  , 
k— ■  1 


(Eqn.  9) 


in  which  j  is 


the  square  root  of  minus  one.  Finally,  ‘the  desired  estimates 
of  the  baseline  vector  and  x„  are  the  values  of  §  and  x  that 


maximize  I C !  =  C*C. 
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With  this  algorithm,  all  of  the  difficulties  of  the 
traditional  algorithm  are  hidden  in  the  problem  of  finding 
the  maximum  of  C*C.  A  priori  information  is  no  less  useful; 
it  serves  to  limit  the  volume  in  4-parameter  space  that 
needs  to  be  searched  to  find  the  maximum.  We  believe  that 
we  can  make  the  search  sufficiently  efficient,  computationally, 
that  this  approach  will  be  feasible.  We  hope  to  present  our 
ideas,  and  some  test  computation  results,  soon. 
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August  12,  1980 


MEMORANDUM 


TO: 

FROM: 

SUBJECT: 

REFERENCE: 


R.  L.  Greenspan 

C.  C.  Counselman  III  and  S.  A.  Gourevitch 

Data  Processing  Algorithms  for  Experiment  I,  Con¬ 
tinued 

Memorandum  from  us  to  you  dated  6/11/80,  "Data 
Processing  Algorithms  for  Experiment  I  ..." 


I.  Introduction 

When  we  left  our  hero,  he  was  estimating  the  unknown  base¬ 
line  vector,  B,  by  analysis  of  the  data  produced  by  the  Mark  III 
processor.  These  data  consisted  of  a  time  series  of  sets  of 
eight  samples  of  the  real-valued  crosscorrelation  function, 
r  ^ ( t . ) ;  here,  the  subscript  i  (=  -4,  -3,  ...,  +3)  is  an  integer 
representing  the  number  of  discrete,  250-nanosecond  "lags"  by 
which  the  signal  from  the  "x"  antenna  is  delayed  in  the  proces¬ 
sor,  relative  to  the  signal  from  the  "y"  antenna;  and  t.  is  the 
i-th  in  the  series  of  times,  uniformly  spaced  with  At  =  t._^  -  t^ 
=  2  seconds,  at  which  the  processor  puts  out  data.  We  had  de¬ 
rived  theoretically  the  expected  form,  ?t(t^),  of  the  correla¬ 
tions  as  a  function  of  a  trial  value,  £,  of  the  baseline  vector, 
and  a  trial  value,  x  ,  of  the  instrumental  group  delay.  The 
algorithm  proposed  for  estimation  of  _B  was  to  maximize  the  mag¬ 
nitude  of  the  correlation,  C,  between  the  observed  correlations, 

r  ,  and  the  theoretical  correlations,  t  .  (Note  that,  whereas 
£  £ 

the  observed  values,  r  ,  are  real-valued,  the  theoretical  values. 


29 


-2- 


r  ,  are  complex;  thus  C  is  complex.)  The  magnitude  |C|  is  a 
fL 

function  of  four  variables:  8,  9»  2,  and  ?  f  of  which  the  first 
three  are  the  Cartesian  components  of  &,  and  the  fourth,  ?  .  is  a 
"nuisance"  parameter.  Fortunately,  |C|  is  a  relatively  weak 
function  of  tq;  thus,  the  essential  computational  problem  is  to 
find  the  maximum  of  a  function  of  only  (!)  three  variables,  x,  5", 
and  z . 

The  purpose  of  the  present  memorandum  is  to  document  our 
further  investigation  of  this  algorithm  for  baseline  estimation. 
For  simplicity,  the  description  will  be  chronological. 

II.  Ambiguity  Mapping 

In  general,  |C|  does  not  have  a  unique  maximum.  In  addition 
to  the  principal,  or  global  maximum  (the  highest)  which,  hope¬ 
fully,  corresponds  to  the  desired  estimate  of  j3,  there  are  many 
secondary,  local  maxima  (with  lesser  values)  that  correspond  to 
possible  spurious  estimates.  We  refer  to  all  the  maxima  as 
"ambiguities,"  and  to  the  process  of  identifying  the  principal 
maximum  as  "resolving  the  ambiguity"  in  the  estimate  of  the 
baseline  vector.  We  call  |C|  the  "ambiguity  function;"  a  map  of 
|  C  |  as  a  function  of  IS  is  called  an  "ambiguity  map." 

To  be  confident  of  finding  the  principal  maximum,  and  cor¬ 
rectly  resolving  the  ambiguity,  one  must  search  a  volume  of  IS 
space  that  encompasses  the  position  of  this  maximum  --  the  de¬ 
sired  estimate  of  B.  This  volume  must  also  be  sampled  at  "grid" 
points  spaced  finely  enough  to  ensure  that  the  desired  estimate 
is  not  overlooked.  Finally,  one  must  ensure  by  the  design  of  the 


30 


-3- 


experiment  that  even  in  the  presence  of  noise  and  other  errors, 
random  and  systematic,  the  height  of  the  maximum  associated  with 
the  desired  estimate  will  exceed  that  of  the  highest  other  max¬ 
imum  in  the  volume  searched. 

In  our  experiment,  additive  white  noise  is  unlikely  to  be  as 
important  as  systematic  errors  due  to  multipath  interference  and 
local-oscillator  phase  drifts.  We  will  discuss  these  error 
sources  in  some  detail  later.  For  the  moment,  we  merely  observe 
that  in  practice  we  will  probably  need  to  search  a  volume  of  the 
order  of  a  few  cubic  meters,  set  by  the  uncertainty  of  available 
a  priori  information  on  the  baseline  vector;  and  the  grid  spacing 
will  be  between  about  10  and  30  centimeters.  The  appropriate 
grid  spacing  is  set  by  the  LI  wavelength,  19  centimeters,  which 
determines  the  characteristic  width  of  the  maxima  of  the  ambigu¬ 
ity  function. 

How  easy  it  will  be  to  distinguish  the  correct  ambiguity  in 
practice  will  depend  not  only  upon  the  magnitudes  of  the  errors 
present,  but  also  upon  the  schedule  of  observations.  The  length 
of  time  spanned  by  the  observations  is  particularly  critical. 

How  long  is  enough?  To  attempt  to  answer  this  question  in  ad¬ 
vance  of  the  performance  of  actual  experiments,  we  have  begun  to 
simulate  experiemnts  in  the  computer.  Figure  1  shows  a  map  of 
| C  j  for  computer-simulated  VLBI  observations  of  the  existing  six 
GPS  satellites  with  a  pair  of  MITES  antennas  separated  by  30 
meters  along  a  horizontal,  East-West  line,  located  near  72°W. 
longitude,  41°M.  latitude.  The  time  span  of  the  observations  was 
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Figure  1. 


Ambiguity  map  for 


1-hour  experiment. 
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Figure  2.  Ambiguity  map  for  2 -hour  experiment 


1  hour. 


This  map  is  of  a  particular  vertical  plane  that  contains 
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the  true  value  of  the  baseline  vector  at  the  origin.  No  signi¬ 
ficant  amount  of  additive  noise,  oscillator  phase  drift,  multi- 
path  interference,  or  other  error  corrupted  the  simulated  obser¬ 
vations,  which  were  of  only  the  Ll  (1575.42  MHz)  signals.  Figure 
2  is  like  Figure  1,  except  the  period  of  observation  is  two  hours 
instead  of  one.  The  normalized  heights  of  the  secondary  maxima 
are  significantly  reduced  by  the  doubling  of  the  period  of  obser¬ 
vation. 

III.  Broken  Coherence 

In  the  simulations  that  yielded  the  ambiguity  maps  shown  in 
Figures  1  and  2,  it  was  assumed  that  the  local  oscillator  phase, 
although  unknown,  was  constant.  We  have  also  simulated  one-  and 
two-hour  experiments  in  which  we  assumed  the  oscillator  phase  had 
one  constant  unknown  value  for  80 0  seconds,  another  statistically 
independent  constant  unknown  value  for  the  next  800  seconds, 
another  independent  value  for  the  next,  and  so  on.  Operation¬ 
ally,  this  meant  that  we  summed  the  complex  correlation,  C,  sepa¬ 
rately  for  each  800-second  sub-interval  of  observations.  The 
magnitudes  of  the  complex  sub-sums  were  then  summed  together  to 
yield  the  value  plotted  in  the  ambiguity  map.  Thus,  the  coher¬ 
ence  of  the  summation  was  broken  every  800  seconds.  Such  a 
procedure  would  be  used  to  process  the  observations  from  a  real 
experiment  if  it  were  believed  that  the  oscillators  remained 
"coherent"  --  that  is,  their  phases  remained  constant  within 


approximately  1  radian  --  for  about  800  seconds,  but  that  they 
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could  not  be  relied  upon  to  remain  coherent  for  a  longer  time. 
(Note  that  a  phase  drift  rate  of  1  radian  per  800  seconds  at 
1575.42  MHz  corresponds  to  a  fractional  offset  in  frequency  of 
about  1.26x10  ^ . )  The  effect  of  this  "broken  coherence"  in  the 
summation  is  illustrated  in  Figures  3  and  4  for  the  1-hour  ana 
the  2-hour  experiments,  respectively.  Mainly,  the  general  level 
of  background  clutter  seems  to  have  been  raised,  and  now  there 
are  many  secondary  peaks  with  heights  about  two-thirds  of  the 
principal-maximum  height.  It  should  also  be  noted  that  the  ratio 
of  the  principal-maximum  height  to  the  peak  background  height  is 
not  much  greater  for  the  2-hour  experiment  than  for  the  1-hour 
one . 

Finally,  in  Figure  5  we  again  display  the  ambiguity  function 
for  the  2-hour  experiment  in  which  the  coherence  of  the  summation 
was  broken  every  800  seconds,  but  in  this  map  the  scale  has  been 
reduced  (the  grid  spacing  has  been  increased)  by  a  factor  of  3, 
relative  to  Figure  4.  Our  purposes  in  showing  this  reduced-scale 
map  are  two:  first,  to  demonstrate  that  no  important  maxima  were 
lurking  just  offstage  in  Figure  4;  and  second,  to  indicate  the 
feasibility  of  using  a  two-step,  coarse  and  fine,  search  pro¬ 
cedure  to  locate  the  global  maximum  in  a  very  large  volume.  A 
volume  of,  say,  200  cubic  meters  might  be  searched  with  a  coarse 
grid  of  30-centimeter  (perhaps  even  larger)  spacing  in  order  to 
define  one  or  more  smaller  regions  of  volume  ^1  cubic  meter,  in 
which  the  ambiguity  function  had  sign i f i can tl y  higher-than- 
ave rage  value.  These  smaller  regions  might  then  be  examined  with 
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Figure  3.  Ambiguity  map  for  1-hour  experiment  with  broken 
coherence.  Circled  points  have  values  above  63% 
of  global  maximum. 
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Figure  4.  Ambiguity  map  for  2 -hour  experiment  with  broken 

coherence.  As  in  Figure  2,  apart  from  the  origin 
there  is  no  grid  point  at  which  the  ambiguity 
function  has  more  than  63%  of  the  global  maximum 
value. 


VI  vl  VI  vl  VI  Ui 


WEST 


6  5  5  6  3 


6  7  6  7 

6  8  6  7  7>V£_ 
6  7  7  ft—6— a— 6  • 
5  4  7  7  6  7  8 

5  6  4  4  5  5  6 

5  4  5  4  5  5  4 

5  4  5  5  3  6  5 

5  5  5  5  4  4  4 

3  4  4  4 

3  3  3  3 

3  3  3  3 

3  3  3  1 

2  3  4  2 

4  3  3  3 

4  3  4  4 


Figure  5. 


7 

9 

8 

8 

81 

¥ 

8 

8 

3 

7| 

[ T 

A 

8 

1a 

Al 

7. 

J 

c 

<L 

B 

B 

a 

"5) 

A 

8 

8 

A 

7 

9 

0 

Vl 

Jj 

0 

\* 

A 

s 

!  J  «  M  -7  J  -9  J  C.  C. 

B  (S  A  Ajs  5  5  3  5  3  3  2 
B/.9  a(T  7S543333 
($,  &)  A^|  7  5  6  3  4  2  2  3 
[S/A  A  C !  9  4743433 
3  C  C  3  7  6  7  3  6  3  4  5 


C  F|9  4 


7  3  5  3  5  4 
7  5  4  6  3  5 


4  3  5  3  3  4 
4  4  5  4  4  4 


l)0  5AEAG  9  7776S33435 
A  S  B  H  E  G  f  ^9/f- B)7  9  6  6  7  7  6  4 
^..5,  C-a-s(Y)s-e-.7.C..4\9-9-5  6-9-8-94- 
6(5  D®F  H  E  H  C  B  B  A*a)  6  6  7  7  7 
6  7\a_A  G  4  F(§)B  0  ®A/9  S  S  6  7  5 
6  7  6  33  D  GFA0C\97389@96 
4  5  9  7|C  B  5  A  F  B  3fi  3(A\o  Y  8  Cl 


6  7  6  8  ]  0  GFA0C\9738 

4  5  9  7jC  B  G  A  F  B  t)6  80 

3  6  3 (A  p  3  A  A  C  Art  9  8  5 

3  5  5  7 " t\0  C  ©  A  8  A)8  6  8 


8  5\05 
6  8  7  6 

6  7  6  6 

7  5  6  6 


2  7  7  7  j)S  C  C  3(8  8  6  76666 

4  4  7  Sfc&Z/'gS*  1)8  7  5  6  6  5  5 

3  6  6  7  A  (Y®  9/3(9  8  6  7  6  6  5  4 

4  4  7  8  A  AV? (X  C|9  9  8  6  7  5  5  5 

3  4  7  aU  A  (sjc  AU  9  8  8  4  6  4  6 

2  5  7  7  4  3^A^^ 7l  o| 6  7  6  7  6  6  4 


LOCAL 

VERTICAL 


’  I  r  -  i  -  i  -  I 

|^—3  meters — J 

Like  Figure  4,  except  the  grid  spacing  is  30 
centimeters.  Only  the  point  at  the  origin  ha: 
value  above  63%  of  the  global  maximum,  which 
is  at  the  origin.  All  grid  points  with  value 
above  23.5%  of  the  global  maximum  are  enclose* 
by  a  solid  contour  in  this  Figure. 
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10-cent imeter  grid  spacing  in  order  to  isolate  the  global  maxi¬ 
mum.  The  final  "peaking  up"  on  this  maximum  would  be  done  by  a 
standard  function-maximizing  routine.  We  have  used  the  routine 
STEPIT  (copyright  1965  by  J.  P.  Chandler,  Physics  Department, 
Indiana  University,  Bloomington)  with  satisfactory  results  for 
these  simulated  experiments. 

IV.  Change  of  Observable:  Complex  Cross-Power 
A  problem  with  ambiguity  mapping,  for  either  a  simulated  or 
a  real  experiment,  is  that  it  can  consume  a  large  amount  of 
computer  time.  For  example,  consider  that  one  cubic  meter  con¬ 
tains  10^  points  of  a  10-centimeter  grid;  to  compute  the  corre¬ 
lation  C  for  one  grid  point  requires  a  summation  over  a  long 
time- series  of  observations  --  perhaps  10^  points  in  time;  and 
for  each  point  in  time,  the  product  r^r^  must  be  summed  over 
several  values  of  i .  Each  of  the  several  million  spectral  den¬ 
sity  function  evaluations  thus  indicated  may,  in  turn,  involve  a 
fairly  large  number  of  elementary  arithmetic  and  logical  opera¬ 
tions.  There  is,  therefore,  some  incentive  to  improve  the  effi¬ 
ciency  of  the  algorithm.  One  improvement  that  we  have  investi¬ 
gated  and  recommend  adopting  is  to  pre-process  the  raw  correla 

tion  data,  the  set  of  eiaht  or  more  r  's,  at  each  point  in  time, 

2, 

in  order  to  obtain  a  single  complex  number  that  embodies  most  of 
the  information  that  the  full  set  of  r  's  contained.  This  com- 

l 

plex  datum,  which  we  denote  by  the  symbol  5,  is  the  cross-power 
spectral  density  function  evaluated  at  the  GPS  carrier  frequency; 
it  is  computed  from  the  r, 's  by  the  Fourier  transform  relation 


?■  7 


( Eqn .  1) 


S  =  ri  *  exp(  » 

where  j  is  the  square  root  of  minus  one.  Thus,  we  compress  eight 
or  more  real-valued  data  into  one  complex  datum  at  each  point  in 
time.  This  reduction  of  the  quantity  of  data  yields  some  saving 
in  computer  time.  However,  an  even  greater  saving  results  from 
the  fact  that  the  theoretical  value  of  S  can  be  computed  more 
easily  than  the  ?  can  be. 

Theoretically,  S  is  related  to  the  k-th  satellite's  received 
power  Ak  and  interferometric  phase  delay  t^,  and  to  the  local- 
oscillator  phase  difference  <t>LQ,  by: 

S  =  (1/2)  ^  Ak*exp(  jwcTk—  j*L0)  *  5  A  (Tk“*6) 

+  (1/2)  ^  Ak*exp(«.j<i)cTk+j<|>L0)  •  (-1)  (xk-i,6)  (Eqn.  2) 


[For  comparison,  refer  to  Eqn.  (1)  of  the  6/11/80  memo;  note  that 

^c’“LO^*4  =  Here,  A(x)  is  the  triangle  function  defined 

in  Equation  (2)  of  the  6/11/80  memorandum.  From  Equation  (2) 

above,  it  may  not  be  obvious  that  the  theoretical  value  of  the 

new  observable,  S,  is  any  easier  to  compute  than  that  of  the  old 

observable,  r^ .  However,  a  good  approximation  to  Equation  (2) 

above  is  given  by  the  simple  relation 

s> 

S  =  2  Ak’  exp(  Du,cTk";)  4,LO)  *  (Eqn*  3) 


That  is. 


^ jr •  4  (  t  l  6  ) 


=  4 


(Eqn.  4) 


and 


?r  (-1)1-A(t k-l«)  = 


0 


(Eqn.  5) 


for  all  k.  These  approximations  are  good  as  long  as  the  range  of 
values  of  i  over  which  the  sums  extend  includes  all  significantly 
nonzero  values  of  the  triangle  function,  Mt^-45).  The  accur¬ 
acies  of  these  approximations  are  discussed  in  detail  in  Appendix 
A. 

The  baseline-estimation  algorithm  set  forth  in  the  6/11/80 
memorandum  is  now  recast  as  follows. 

1.  First,  condense  the  set  of  observed  values  r  (t.)  at 
each  point  in  time  t^  to  a  single  complex  number  S(t^). 

2.  Compute  the  complex  correlation  C  between  the  "observed" 
function  of  time  S(t-),  and  a  trial,  theoretical,  function  ^(tj) 
by 

^  \  Ak.exp(  ju)ctk)  ;  ( Eqn .  6) 

C  =  ^  S(t.)t*(t.)  .  (Eqn.  7) 

The  summation  over  i  extends  over  the  entire  time  span  of  the 
observations  if  the  local  oscillator  phase  is  assumed  to  remain 
constant  for  the  entire  span.  Otherwise,  C  may  be  sub-summed 
over  sub-spans  as  described  in  Section  III. 

3.  Take  the  magnitude  of  C,  or  the  sum  of  the  magnitudes  of 
the  complex  sub-sums,  and  find  the  baseline  vector  that  maximizes 
this  magnitude. 


This  recast  algorithm  is  faster  than  the  one  originally 
proposed  by  a  factor  of  about  4.  The  improvement  in  speed  comes 
partly  from  the  elimination  of  any  summation  over  i  in  Equation 


(7),  and  partly  from  the  elimination  of  the  triangle  function  in 
Equation  (6).  To  obtain  the  improved  speed,  we  have  committed  a 
conscious  error  by  substituting  the  simpler  form  (3)  [which  leads 
to  (6)]  for  the  more  exact  form  (2). 

What  is  the  effect  of  this  error  on  the  ambiguity  map,  and 
on  the  estimate  of  the  baseline  vector?  For  the  simulated  exper¬ 
iments  described  in  Section  II,  we  have  re-computed  the  ambiguity 
function  --  the  magnitude  of  C  --  using  Equation  (7).  In  these 
experiments  the  range  of  values  of  i  for  which  the  r^  were  ob¬ 
tained  was  i  =  -4  to  +3.  We  find  that  the  largest  magnitude  of 
the  difference  between  the  ambiguity  function  computed  "exactly," 
according  to  the  6/11/80  memorandum,  and  that  computed  according 
to  Equation  (7)  above,  is  about  3%  of  the  principal  maximum 
value.  The  position  of  the  principal  maximum  is  not  changed  by 
more  than  a  few  millimeters,  which  is  the  level  of  uncertainty 
due  to  our  use  of  single-precision  arithmetic  in  the  computa¬ 
tions.  We  conclude  that  the  use  of  the  faster  algorithm  does  not 
significantly  increase  the  probability  of  mistaking  a  secondary 
for  the  principal  maximum,  or  degrade  the  accuracy  of  the  base¬ 
line  estimate  given  by  the  position  of  the  principal  maximum.  If 
the  slight  loss  of  position  accuracy  should  be  deemed  important, 
one  could  always  refine  the  baseline  estimate  by  reversion  to  the 


original  algorithm,  after  having  obtained  a  preliminary  estimate 
by  means  of  the  faster  algorithm.  Note  also  (see  Appendix  A) 
that  the  accuracy  of  the  faster  algorithm  would  have  been  better, 
had  the  r .  data  been  available  for  a  wider  range  of  l . 


The  triangle  function,  a(x),  is  defined  by 


f  1  -  1.023  MHz*  j  x |  ) 

A ( x )  =  J  or  ^  ,  whichever  is  greater. 

This  function  approximates  the  autocorrelation  function  of  the  C/A 
modulation  of  the  GPS  signal.  For  convenience,  we  also  define  £  + 
and  Z~  by 


max 

z  +  =  ^/\  A  (  t  —  2.  6  )  , 


i  =i 


min 


and 


lmax 

(-1) t*A{T-t6) , 

where  6  =  250  nanoseconds.  In  the  text  it  is  claimed  that  Z+  =  4 
and  z~  =  0  as  long  as  A  (  t  — <5 )  =  0  for  and  for  £>£max*  Th6 

accuracies  of  these  two  approximations  may  be  judged  by  inspection 
of  Figures  A.l  and  A. 2.  If  the  stated  condition  on  the  range  of  i 
is  satisfied,  then  r+  =  3.91  and  z  =0  (exactly).  As  long  as  z 
vanishes  exactly,  then  there  is  no  loss  of  geodetic  accuracy  if 
has  any  constant  nonzero  value.  Thus,  the  difference  between  3.91 
and  4  is  of  no  significance. 
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I.  Introduction 

The  geodetic  interferometry  experiments  that  we  are  planning 
to  do  with  antennas  of  the  MITES  design  (crossed  horizontal 
dipoles,  fed  in  phase  quadrature,  located  3A/8  above  and  oriented 
parallel  to  a  metallic  "ground"  plane)  will  have  results  affected 
to  some  degree  by  interference  from  signals  reflected  from  nearby 
surfaces  such  as  those  of  the  earth,  building  walls,  and  fences. 
We  are  attempting  to  predict  the  effects  of  this  interference  by 
using  the  computer  to  simulate  observations  with  various 
geometries,  time  spans,  etc.  The  mathematical  formulas 
underlying  our  computer  code  are  derived  in  this  memorandum. 
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II.  Vector,  Matrix,  and  Operator  Notation 

We  use  the  following  notation: 

Scalar  quantities  are  denoted  by  either  lower-  or  upper-case 
symbols,  for  example  p  or  E. 

A  vector,  which  has  a  magnitude  and  a  direction,  is  denoted 
by  a  lower-  or  an  upper-case  symbol  with  an  arrow  above  it,  for 
example  p  or  E.  The  magnitude  of  any  vector  is  represented  by 
the  same  symbol  without  the  arrow.  The  unit  vector  in  the  same 
direction  is  represented  by  the  same  symbol  with  a  "hat"  instead 
of  an  arrow.  Thus,  p  =  p£>.  Two  other  symbols  for  a  vector,  <p] 
and  [p>,  will  be  introduced  below  in  order  to  facilitate  the 
writing  of  outer  or  tensor  products.  In  this  memorandum  just  one 
coordinate  system  is  used  for  the  expression  of  vectors  in 
component  form,  and  it  is  sufficient  for  us  to  define  <p]  simply 
as  the  row,  and  [p>  as  the  column,  matrix  of  the  components  of 
the  vector  p.  A  fuller  explanation  of  the  meanings  of  these 
symbols  is  given  in  the  Appendix. 

Scalars  and  vectors  may  be  complex .  The  complex  conjugate 
is  formed  by  "starring".  For  example,  p*  is  the  conjugate  of  p. 

A  right-handed  Cartesian  coordinate  system  will  be  used 
whose  x  and  y  axes  are  parallel  to  the  crossed  dipoles  of  a  MITES 
antenna;  the  positive  z  axis  extends  upward.  Unit  vectors  along 
the  positive  axis  directions  are  x,£,  and  z. 

The  3x1  column  matrix  of  the  x,  y,  and  z  coordinate 
components  (each  a  scalar)  of  a  vector  v  is  denoted  by  [v>.  That 


-3- 


The  1x3  row  matrix  of  the  Cartesian  components  of  v  is 
denoted  by  <v]  . 

The  "inner",  or  scalar  product  u’v  of  two  vectors  u  and  v  is 
therefore  given  by  the  matrix  product 
U*  V  =  <u]  [ v> 

which  is  also  written  as  <u|v>. 

The  "outer",  or  tensor  product  of  u  and  v  is  represented  by 
[u><v]  . 

The  parallel  projection  operator  Pu  is  a  square  matrix 
defined  for  any  nonzero  vector  u  by 
Pu  =  [ u>  <u] . 

Pre-multiplication  by  Pu  of  the  column  matrix  [v>  yields  the 

column  matrix  of  the  components  of  a  new  vector  which  is  parallel 

to  u  and  whose  magnitude  is  equal  to  the  magnitude  of  the 

projection  of  v  along  u.  Thus,  Py  "projects"  a  vector  along  u. 

It  is  obvious  from  this  description,  and  may  easily  be  shown  by 

direct  multiplication  (since  <u|u>  =  1),  that  PuPu  =  P*  =  Py. 

That  is,  the  parallel  projection  of  the  parallel  projection  is 

the  same  as  the  parallel  projection. 

The  perpendicular  projection  operator  Q(J  is  defined  by 

Q  =  I  -  P  , 
vu  u 

where  I  is  the  3x3  identity  matrix.  We  call  this  operator  a 
perpendicular  projection  operator  because  it  subtracts  away  the 
parallel  projection  of  the  vector  upon  which  it  operates,  leaving 
that  part  of  the  vector  which  is  perpendicular  to  u.  Note  that 


0  0  =  Q2  =  Q  . 

*;uvu  vu  u 

The  reflection  operator  Ru  is  defined  by 


Ru  "  1  "  2V 
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For  simplicity  in  the  algebraic  manipulations  that  follow,  we 


assume  that  p^  =  1.  Note,  however,  that  with  p  =  1,  the 
magnitude  of  the  complex  "length"  of  the  vector  p  is  not  equal  to 
unity.  The  square  of  the  magnitude,  given  by  <p*|p>,  equals  2. 


IV.  Effect  of  Reflecting  Plane 
We  specify  the  location  and  the  orientation  of  a  reflecting 
plane  by  a  "mirror"  vector  m  which  is  just  the  position  vector, 
with  respect  to  the  origin  of  coordinates,  of  the  point  in  the 


plane  that  is  nearest  the  origin.  Thus,  the  mirror  vector  is 
perpendicular  to  the  plane  and  points  toward  it.  In  terms  of 
this  vector  it  is  easy  to  compute  the  electric  field  vector  of 
a  reflected  wave  in  terms  of  the  field  of  an  incident  wave.  We 
know  that: 


1.  The  reflected  wave  vector  is  related  to  the  incident 

wave  vector  k ^  through  the  reflection  operation: 

|  k  >  =  R  I  k  .  > 

'  r  m  i 


and 


'V  -  'V- 

The  reflected  wave  that  arrives  at  a  far-distant  position 

r  traveling  in  the  direction  of  r  was  originally  radiated 

in  the  direction  of  the  unit  vector  k.  whose  components 

are  given  by  R  lr>.  (Here  we  assume  r>>m,  and  ianore 
m 

parallax.)  Therefore  in  Equation  1 1 1  —  1 , 

Qk  =  I-Pk  =  I-lkXk|  =  I  -  Rm  lr><r|Rm. 

Since  RJ  =  I,  we  may  write 
m  J 
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2.  At  the  reflecting  plane  the  incident  and  the  reflectec 
fields  are  related  through  a  boundary  condition.  For 
simplicity  in  this  memorandum  we  assume  that  reflecting 
planes  are  perfectly  conducting  so  that  the  sum  of  the 
in-plane  components  of  the  electric  fields  (direct  plus 
reflected)  is  zero.  Thus,  the  field  is  reflected  with  a 
reversal  of  the  sign  of  the  in-plane  component  and  no 
change  of  the  normal  component.  At  the  plane,  therefore, 
we  have 

'V  ■ 

Combining  this  result  with  the  result  for  Qk  derived  in 
paragraph  1  above,  we  obtain 

I E r >  =  -Rm  ‘(scalar  factors)*  RmQrRmlF» 

=  -  (scalar  f actor s) • QrRm I P> , 
where  the  "scalar  factors"  depend  only  on  the  path 
length  . 

3.  The  length  of  the  reflected  ray  path  exceeds  that  of  the 
direct  ray  path  by  2k^*m  =  -2k^*m.  (Again,  we  are 
ignoring  parallax.)  This  addition  to  the  path  length  may 
be  ignored  in  the  r  2  factor  because  r>>m,  but  must  be 
incorporated  in  the  complex  exponential  phase  factor 
because  generally  m  is  not  negligible  compared  to  the 
wavelength.  The  added  path  length  is  accounted  for  by 
substituting  (r-2m)  for  r  in  the  argument  of  the  complex 
exponential. 

Now  combining  the  results  cf  paragraphs  1,  2,  and  3,  we  have 
|Er>  =  -k  2  r  ~  1  exp( i<k I  ( r  —  2  if. )  >  )  QkRmlp>  , 


48 


in  which  £  =  kr.  Thus,  the  expression  for  the  reflected  field  at 
a  distance  r>>m  is  the  same  as  the  original  expression  derived  in 
Section  III  but  with  two  changes: 

1.  The  dipole  moment  | p>  is  replaced  with  its  negative 

reflection,  -R  |  p>  . 

'  m  r 

2.  The  vector  r  representing  the  position  of  the 

observer  with  respect  to  the  antenna  at  the  origin  is 
+  •> 

replaced  by  (r-2m),  which  is  the  observer's  position 
with  respect  to  the  image  of  the  antenna  in  the 
reflecting  plane. 


V.  Effect  of  Horizontal  "Ground”  Plane 
In  this  Section  we  apply  the  result  derived  in  Section  IV  to 
compute  the  effect  of  the  MITES  antenna's  ‘'ground"  plane,  which 
is  described  by  the  mirror  vector 

r  o  i 


lg> 


o 

-g 


where  g  =  3X/8.  At  a  distance  r>>g  from  the  origin,  the  electric 

field  of  the  wave  reflected  by  the  ground  plane  is 

IE  >  =  -k2r_1exp{ i<k ! (r-2g) >)  Q.  R  I p> ,  (Eqn.  V-l) 

r  k  cj 

where  k  =  kr. 

■4- 

But  because  for  the  MITES  antenna  p  is  perpendicular  to  g, 

=  lP>. 

Thus,  in  this  special  case  (only.1),  the  operator  RQ  may  be 
deleted  from  Equation  (V-l).  With  this  deletion  the  expression 
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for  the  reflected  field  becomes  very  similar  to  the  one  for  the 
direct  field,  and  the  two  are  easily  combined  to  obtain  the  total 
fieid,  Etot: 

IE  t>  =  k2 t~x exp( i<£ I r> ) (1 -exp( - 2i <£ I g> ) Qk i P> . 

This  expression,  in  turn,  is  simplified  by  moving  the  origin  of 
coordinates  to  the  point  in  the  ground  plane  directly  below  the 
center  of  the  crossed  dipoles  --  that  is,  to  the  point  with  "old" 
coordinates  |g>.  In  "new”  coordinates,  the  total  field  radiated 
by  a  MITES  antenna  is 

IE  >  =  k2 r~ Xexp( i< k | r > ) ( exp( i< k | g> ) -exp( -i < k ] g> ) ) Q.| p> 

=  2ik2r  ^exp { i<k | r > ) s i n ( <k I g> ) Qk I p> ,  (Eqn.  V— 1 ) 

where  again  we  have  neglected  the  change  in  the  r  2  factor. 


VI.  Signal  Received  from  a  Transmitter 
The  complex  amplitude  of  the  signal  received  by  an  antenna 
at  a  distant  position  r  is  proportional  to  the  inner  or  "dot" 
product  of  the  complex  conjugate  of  the  receiving  antenna's 
effective  electric  dipole  moment  vector  with  the  radiated 
electric  field  vector  evaluated  at  r.  Denoting  the  receiving 
antenna's  electric  dipole  moment  by  q,  we  have 

A  =  <q*|£  >  (Eqn.  VI-1) 

d 

where  A  is  proportional  to  the  received  signal  amplitude.  By 
"effective"  dipole  moment  here  we  mean  the  dipole  moment  of  an 
infinitesimal  antenna  having  the  same  polarization  as  the  actual 
antenna,  multiplied  by  the  "gain"  of  the  actual  antenna.  The 
gain  factor  can  account,  for  example,  for  the  effect  of 
incorporating  a  reflector  or  multiple  array  elements  into  one 


"antenna"  structure.  By  "signal"  here  we  mean  the  current  or  the 
voltage  delivered  to  a  load.'  The  constant  of  proportional  i  ty 
between  A  and  the  received  signal  amplitude  contains  not  only  the 
antenna  gain  factor  but  also  various  impedances,  etc.  Assuming 
that  it  _i_s  constant,  we  are  not  concerned  with  it  and  will  not 
discuss  it  further. 

The  MITES  antennas  will  be  used  to  receive  signals 
transmitted  by  the  GPS  satellites,  not  to  radiate  signals  for  the 
satellites  to  receive.  But  for  expository  purposes  in  this 
memorandum  we  have  been  calculating  the  field  radiated  by  a  MITES 
antenna  as  if  J_t  were  transmitting.  We  continue  in  this  vein, 
and  now  calculate  the  amplitude  of  the  signal  that  a  GPS 
satellite  antenna  would  receive  from  a  MITES  antenna.  By 
reciprocity,  we  know  that  the  result  for  the  received  amplitude 
will  be  the  same  as  we  would  obtain  if  we  treated  the  GPS  antenna 
as  the  radiator  and  the  MITES  antenna  as  the  receiver. 

A  GPS  satellite  antenna  has  circular  polarization  with  q 
perpendicular  to  r  and  £.  (GPS  satellites  are  always  oriented  so 
that  their  antennas  face  earthward.)  Thus,  q  can  be  expressed  as 
a  simple  rotation  of  p: 

I q>  =  D | p>  ,  ( Eqn  .  VI-2) 

where  D  is  a  matrix  rotation  operator  that  depends  on  the  azimuth 

4>  and  the  zenith  angle  6  of  the  satellite  relative  to  the  MITES 

antenna.  The  matrix  D  can  be  written  as  the  product  of  two 

matrices  D  and  D  ,  the  first  beino  a  function  of  the  azimuth 
0  6 

alone,  and  the  second  a  function  of  only  the  zenith  angle: 


represents  a  rotation,  about  the  y  axis,  of  the  z  axis  toward  the 
x  axis,  and 

cos4>  -sin$  0 

D  =  sin$  cos4>  0 

$ 

0  0  1 

represents  a  rotation  about  the  (original)  z  axis,  of  the  x  axis 
toward  the  y  axis.  Thus, 

cos6cos<t>  -sin<t>  sin0cos<t> 

D  =  cosSsimt’  cos  4>  sin0sin<t> 

-sine  0  cose 

(Note  that  we  are  not  accounting  for  a  possible  third  rotation, 
about  an  axis  parallel  to  the  wave-propagation  direction  r.  In 
other  words,  we  are  ignoring  the  "third"  Euler  angle,  \|>.  Such  a 
rotation  would  only  change  the  phase  of  A,  by  the  same  angle  as 
the  rotation,  by  virtue  of  the  circularity  of  the  polarization  of 

,  i  \h 

the  GPS  antenna.  That  is,  D  =  e  .) 

From  the  definition  of  D  it  follows  that 

I  k>  =  D  | z> 
and  . 

Qk  =  DQZD  .  (Eqn .  VI-3) 

Combining  Equations  (V-l),  (VI-1),  (VI-2),  and  (VI-3),  and 

2-1 

omitting  the  constant  factor  (2ik  r  ) ,  we  obtain 

_  i  ikr 

A  =  <p*IQ2D  | p>  sin(k*g)e 

-*■  _  l  4^  jkr 

=  <p* ID  lp>  sin(k’g)e 

•4  -4  4  4  I  r 

=  < p 1 D | p*>  sin  (k*  g) e 


(Eqn.  VI-4) 
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By  direct  substitution  and  multiplication  we  find 

< p | D | p* >  =  e 1  *  (1  +  cose).  (Eqn.  VI-5) 

The  received  power  is  proportional  to 

A*A  =  (1+cose)2  sin2(<£|g>) 

2  2 

=  (1+cose)  sin  ( (3»/4) cose ) .  (Eqn.  VI-6) 

Thus  --  except  for  a  constant  factor  --  we  have  derived  the 
directive  power  gain  of  the  MITES  antenna  as  a  function  of  zenith 
angle,  for  circular  polarization.  The  power  gain,  as  expected, 
varies  with  zenith  angle  but  is  independent  of  azimuth.  Equation 
(VI-5)  shows,  however,  that  the  received  signal  phase  varies  with 
azimuth  --  in  fact,  is  just  equal  to  the  azimuth  --  whereas  it  is 
independent  of  zenith  angle. 

The  constant  factor  by  which  the  right  side  of  (VI-6)  must 
be  multiplied  in  order  to  obtain  the  power 'gain  of  the  MITES 
antenna  with  respect  to  the  standard  (fictitious),  isotropic, 
circularly  polarized,  reference  antenna  is  determined  by  the 
condition  that  the  integral  of  the  gain  over  the  visible 
hemisphere  (0^^/2;  b^<2*  ;  note  that  the  gain  is  zero  for  the 
other  hemisphere,  equals  4n.  We  find  that  the  numerical 

value  of  the  factor  is  about  1.23.  However,  for  engineering 
purposes,  and  especially  considering  that  we  have  neglected  the 
effects  of  dissipation,  the  finite  sizes  of  the  dipole  elements 
and  the  ground  plane,  etc.,  the  value  of  the  factor  is 
approximated  well  enough  by  unity.  The  function  of  6  that 
appears  on  the  right  side  of  (VI-6)  will  be  denoted  by  G(6): 

G ( e )  =  (1+cose)2  sin2 (  (3^/4) cose ) ,  (Eqn.  VI-7) 

and  we  will  refer  to  G(e)  simply  as  the  "antenna  gain."  Values 
of  this  function  are  given  in  Table  1. 
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Table  1.  Directive  Power  Gain  (ref.  isotropic,  circ.  pol.)  of 
MITES  Antenna  vs.  Zenith  &  Elevation  Angles,  from  Equation 
(VI-7)  . 


Zenith 

Angle ,  e 
(deg) 

Elevation 

Ang  1  e 
(deg) 

Gain,  G 
( d imens ion¬ 
less) 

10  log1Q 

(dbi  c ) 

0 

90 

2.00 

+  3.01 

10 

80 

2.11 

3.24 

20 

70 

2.41 

3.82 

30 

60 

2.77 

4.42 

40 

50 

2.95 

4.70 

50 

40 

2.69 

4.30 

60 


30 


1.92 


+  2.83 
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VII.  Interferometric  Cross-Power  Observable 
The  interferometric,  complex,  cross-power  observable  S  that 
was  defined  in  our  August  12,  1980,  memorandum  (revised  August 
13,  1980)  for  a  single  GPS  satellite,  except  for  constant 
factors,  is  given  by 
S  =  B*A 

where  A,  given  by  Equation  (VI-4),  is  the  complex  amplitude  of 
the  signal  received  from  the  MITES  antenna  located  at  the  origin 
of  coordinates,  and  B  is  the  complex  amplitude  of  the  signal 
received  from  an  identical  MITES  antenna  located  at  the  position 
b.  (We  continue  to  use  the  "new"  definition  of  the  origin  of 
coordinates  that  we  introduced  in  Section  V:  the  point  in  the 
ground  plane  directly  below  the  center  of  the  crossed  dipoles  of 
the  MITES  antenna.)  To  obtain  an  expression  for  the  complex 
amplitude  B  we  simply  take  the  expression  (VI-4)  for  A  and 
substitute  r-S  for  r  in  the  argument  of  the  complex  exponential 
phase  factor  elkr  =  exp(i<k|r>).  We  ignore  the  parallax  of  the 
baseline  —  that  is,  we  ignore  the  difference  between  the 
directions  of  the  wave  vectors  at  the  two  MITES  antenna  locations 
—  and  we  assume  that  the  two  antennas  have  parallel 
orientations,  so  that  their  dipole  moment  vectors  and  hence  their 
azimuth  phase  factors  are  identical.  Thus,  we  obtain 
S  =  B*A  =  G ( e )  exp(i<k|b>) 
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VIII.  Multipath  at  One  Antenna 
The  reflection  of  the  field  by  the  ground  plane  at  each 
MITES  antenna,  which  we  considered  in  Section  V,  is  an  example  of  j 

multipath  propagation,  but  it  is  an  uninteresting,  or  "benign," 
example  so  long  as  the  ground-plane  mirror  vector  g  and  the  ; 

dipole  moment  vector  p  are  the  same  at  each  antenna.  Multipath 


is  interesting,  or  "pathological,"  only  when  it  is  present  at  one 
antenna  in  a  form  that  is  not  duplicated  at  the  other  antenna. 

To  simulate  pathological  multipath  we  assume  that  at  one  antenna 
—  the  one  at  the  origin  of  coordinates  --  in  addition  to  the 
horizontal  ground-plane  reflector  there  is  a  vertical  plane 
"mirror,"  described  by  a  horizontal  mirror  vector,  m,  given  by 


[m>  = 


m  cosy 


m  siny 


By  virtue  of  the  reflection  in  this  mirror  the  complex  amplitude 

A  of  the  signal  received  at  the  satellite  from  the  original  MITES 

antenna  is  augmented  by  A  .  To  derive  an  expression  for  A  we 

m  m 

apply  the  arguments  of  paragraphs  1,  2,  and  3  of  Section  IV  to 
Equation  (VI-4),  and  obtain 


Am  =  -<p* I D-iRm | p>  sin {<£ I Rm I g>)  exp( i<£ I (r-2m) >) , 

( Eqn .  VIII-1) 

where,  as  in  Equation  (VI-4),  we  have  omitted  the  constant  factor 
2  -1 

(2ik  r  ).  Since  the  ground  plane  is  horizontal  and  the  other 
mirror  is  vertical,  <g|m>  =  0  and  <  k  |  | g  >  =  <k|g>.  Thus,  in 

this  special  case,  Equation  (VIII-1)  reduces  to 

Am  =  -<p* I D~ 1 | p>  sin(<k|g>)  exp ( i <k I  ( r-2m)  > ) , 


(Eqn.  VI 1 1- 2 ) 
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Now  the  interferometric  cross-power  observable  is 

S .  .  =  S  +  S 

tot  m 

where 

S  =  S (b)  =  B*A  =  G ( 6 )  exp ( i < k I b> )  (Eqn.  VIII-3) 

is  due  to  the  signal  received  directly  (without  reflection  in  the 
vertical  mirror),  and 

=  B*A^  =  —  { G  {  6  )  -2  M  s i  n'- ( < k  |  g>  )  }  exp ( i < k I ( b-2m) > ) , 

(Eqn.  VI 1 1-4 ) 

in  which 

M  =  <p*  !D|pXp|mXm|D|p*>  .  (Eqn.  VIII-5) 

In  deriving  Equations  (VII 1—4)  and  (VIII-5)  we  substituted 
I  -  2P^  =  I  -  2|mXm|  for  the  operator  in  Equation  VIII-1,  and 
used  <p*  I  D  ^  |  fn>  =  <  rh  I  D  |  p*  >  .  Substitution  of  the  definitions  of 
m,  D,  and  p  into  Equation  VIII-5  yields 

M  =  (1+cose )  -  sin2e  cos(u-<tO  e^w  ^  ^ 

=  {( 1+cose ) 2/2 }{ 1 -[ tan2 ( e/2 )) e2i ^ v  ^  } 

(Eqn.  VII 1-6) 

The  quantity  M  can  also  be  written  in  a  form  that  involves  only 
dot  products: 

M  =  1  +  <k|z>  -  <p*|k>  <k|m>  <fn|p> 

Finally,  substitution  of  (VIII-6)  into  (VIII  —  4)  yields  the  simple 
result 

Sm  =  S ( b-2m) *  F ,  (Eqn.  VIII-7) 

where 

F  =  - [ tan  2 ( e /2 ) )  e2l(u_*}  .  (Eqn.  VIII-6) 

Note  that  the  magnitude  of  F  is  a  function  of  only  the  zenith 
angle ,  e  ,  whereas  the  phase  angle  depends  on  only  the  azimuth  of 
the  satellite  relative  to  the  mirror. 
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The  magnitude  dependence  on  zenith  angle  is  easy  to 
understand  physically:  In  -the  zenith  direction  the  KITES  antenna 
has  circular  polarization  matched  to  the  wave  that  propagates  to 
it  directly  from  the  satellite.  But  the  sense  of  circular 
polarization  of  a  wave  is  reversed  by  reflection  at  a  perfectly 
conducting  surface,  so  the  MITES  antenna  has  zero  response  to  a 
reflected  wave  coming  from  the  zenith.  At  the  zenith,  IF!  =  0. 
Viewed  from  the  horizon,  the  MITES  antenna  has  horizontal,  linear 
polarization  so  that  it  responds  equally  to  waves  with  the  two 
opposite  senses  of  circular  polarization.  At  the  horizon,  t  F I  = 
1.  The  form  of  the  e-dependence  between  the  zenith  and  the 
horizon  may  also  be  derived  directly  by  detailed  consideration  of 
the  antenna  polarization.  Since,  in  most  practical  observing 
situations,  reflections  would  be  more  likely  to  come  from  near 
the  horizon  than  from  near  the  zenith,  we  expect  that  usually  the 
magnitude  of  F  will  be  near  unity. 

The  phase-angle  dependence  on  the  azimuth  of  the  satellite 
relative  to  the  mirror  is  also  easily  understood  physically. 

This  dependence  stems  from  the  e1^  term  in  Equation  (VI-5):  The 
phase  of  the  signal  received  by  a  MITES  antenna  changes  with  the 
azimuth  from  which  the  wave  approaches  the  antenna.  But  this 
azimuth  is  changed  by  2  ( u  —  ♦ )  by  the  reflection  in  the  vertical 
mirror. 

The  negative  sign  of  the  right  side  of  (VIII-8)  is  due  to 
the  reversal  of  the  sign  of  the  horizontal  component  of  the 
electric  field,  that  occurs  upon  reflection  at  the  perfectly 
conducting  surface  of  the  vertical  mirror. 
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IX.  Proaram  Loaic 


In  this  Section  we  outline  the  logic  of  the  computer  program 
that  simulates  observations  in  the  presence  of  multipath. 

At  each  time  for  which  an  observation  is  to  be  made,  the 
program  executes  a  "DO"  loop  over  all  of  the  GPS  satellites. 
Before  this  loop  is  entered,  the  total  value  of  the 
interferometric  observable  is  set  equal  to  zero.  Within  the 
loop,  the  contributions  of  the  individual  satellites  to  the 
observable  are  added  into  this  total.  A  satellite  may  contribute 
not  at  all  (if  it  is  below  the  horizon  or  occulted  by  the 
vertical  mirror),  via  direct  propagation  alone,  or  through 
multipath  propagation,  depending  upon  its  direction  relative  to 
the  observer's  horizon  and  the  vertical-mirror  plane,  and  upon 
the  size,  shape,  and  distance  of  the  mirror. 

For  each  satellite,  the  unit  vector  r  in  the  direction  of 
the  satellite  is  computed  from  the  orbital  elements,  etc.  The 
sign  of  the  dot  product  of  r  and  the  unit  vector  z  in  the 
direction  of  the  observer's  local  vertical  is  tested  to  see 
whethei  the  satellite  is  above  the  observer's  horizon.  If  it  is 
not,  then  no  further  computation  is  done  for  this  satellite. 
Otherwise,  the  dot  product  of  r  with  the  unit  vector  m  in  the 
direction  of  the  mirror  (see  Section  VIII)  is  tested  to  see 
whether  the  satellite  is  occulted  by  the  vertical  plane  mirror. 

We  consider  this  mirror  to  have  a  finite  size;  its  shape  is 
semi-circular,  like  one  end  of  a  quonset  hut.  The  cosine  of  half 
the  angle  subtended  at  the  observing  antenna  by  the  diameter  of 
the  mirror  is  specified  by  a  program  input  parameter  named  QUON. 


If  <rjm>  is  greater  than  QUON,  then  the  satellite  is  occulted  by 
the  mirror  and  no  further  computation  is  done  for  this  satellite. 
Otherwise,  the  value  S(B)  of  this  satellite's  "direct" 
contribution  is  added  to  the  interferometric  observable  (see 
Section  VII).  This  contribution  is  due  to  the  wave  propagated 
directly  to  the  antenna  from  the  satellite.  Now  we  test  <rim> 
again  to  determine  whether,  from  the  position  of  the  antenna,  an 
image  of  the  satellite  appears  within  the  diameter  of  the  mirror. 
If  <r|m>  is  greater  than  (-QUON) ,  then  no  such  reflection  is 
visible  and  no  further  computation  is  done  for  this  satellite. 
Otherwise,  the  value  of  the  interferometric  observable  is 
augmented  by  =  S(B-2m)  •  F,  the  contribution  due  to  the  wave 

reflected  to  the  antenna  by  the  mirror  (see  Section  VIII).  This 
completes  the  computation  for  one  satellite  at  one  time. 

We  conclude  by  correcting  a  minor  mis-statement  that  we  made 
deliberately,  in  order  to  simplify  the  preceding  description: 

What  our  simulation  program  computes  at  each  observing  time  is, 
in  reality,  not  the  complex  cross-power  S,  but  the  set  of  eight 
real-valued  cross-correlations  r  ,  t=-4 ,-3, . . . ,+3 ,  defined  in  our 
June  11,  1980,  memorandum.  For  direct-path  propagation,  we 
compute  rfc  for  the  baseline  B  according  to  Equation  (1)  of  that 
memorandum.  For  the  ref lected-path  propagation,  we  compute  r£ 
for  the  baseline  b-2m,  and  multiply  the  result  by  F,  by  analogy 
with  Equation  (VIII-7)  in  the  present  memorandum.  It  is  easy  to 
see  that  this  analogous  procedure  is  the  correct  one,  because  S 
and  r  are  linearly  related.  (See  our  memorandum  dated  August 
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Appendix  on  Vector  Notation 

To  facilitate  the  writing  of  outer  or  tensor  products  we 
introduce  two  special  symbols,  <u]  and  [v>,  that  represent 
vectors  in  two  different,  but  complementary  spaces.  The  spaces 
are  related  such  that  a  vector  product  is  formed  by  multiplying 
together  one  vector  from  each  space.  Thus  the  inner  or  "dot" 
product  u'v  =  v • u  =  <u!v>  =  <v|u>,  and  the  outer  or  tensor 
product  of  v  and  u  is  written  [vXu].  The  symbols  <u)  and  [v> 
behave  as  row  and  column  matrices,  respectively,  and  indeed,  once 
a  coordinate  system  is  adopted,  the  vector  components  may  be 
written  in  matrix  form  as  in  Section  II  of  this  memorandum.  If 
we  write  [i>  =  [x>,  [y>,  and  [z>  for  i  =  1,  2,  and  3, 

respectively,  then  the  identity  operator  I  can  be  written 
3  -  . 

I  =  i  [  i  >  <  i  ]  ; 
i  =  l 

tv>  =  xfi><i|v>  =  iti>vi 
i  i 


thus , 


<  u  ]  =  I <  u | i >  < i )  =  x  < i 1 u . 

i 

<u|v>  =  <u  I  I  !  v>  =  I  <  u  I  i  ><  i  lv>  =  tu^j  =  u'v. 

i  i 


The  components  of  an  operator  M  are  M—  =  <iiM|j>  ,  so  that 

components  of  the  expression  (v>  =  M[u>  are  found  from 

v.  =  <i|v>  =  x<i|M|j><jlu>  =  x  M .  .  u  .  , 

i  1 1  D 

i  j 

whereas  <v]  =  <u]M  has  components 

v.  =  <  v  I  i  >  =  l<u!jXj|Mi>  =  Xu.r^.. 

a  j  J 

Evidently,  N  is  the  transpose  of  ^ . 
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Memorandum 

To:  R.  L.  Greenspan 

From:  C.  C.  Counselman  III 

» 

Subject:  Comparison  of  GPS  Radio  Interferometry  Determination  of 

a  Triangle  of  Baselines  with  Independent  Determination 
by  Conventional  Survey 

I ntroduction 

On  December  17,  1980,  we  set  a  portable  "Miniature 
Interferometer  Terminal  for  Earth  Surveying"  (MITES)  antenna  atop 
each  of  three  temporary  survey  marks  that  we  had  established  near 
the  Haystack  Observatory  in  Westford,  Massachusetts.  These 
antennas  yielded  radio  interferometric  observations  of  the 
NAVSTAR  Global  Positioning  System  (GPS)  satellites.  We  analyzed 
the  data  from  an  observation  period  of  about  78  minutes,  by  the 
method  decribed  by  C.  C.  Counselman  and  S.  A.  Gourevitch  in  the 
October  1981  issue  of  IEEE  Transactions  on  Geoscience  and  Remote 
Sensing,  to  determine  the  baseline  vectors  between  the  marks.  We 
repeated  the  experiment  on  December  29,  1980,  with  a  different 
antenna  on  each  mark.  On  the  latter  date,  we  also  performed  an 
experiment  to  determine  how  severely  the  baseline  determination 
would  be  affected  by  the  parking  of  a  car  next  to  one  antenna 
during  the  observing  period.  My  Dasher  station  wagon  was  parked 
broadside  to  one  antenna,  with  its  side  about  60  cm  to  the  south¬ 
east.  Observations  made  under  this  condition  were  analyzed  to 
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obtain  a  second  set  of  baseline  de te rm i na t i ons  for  this  day.  The 
results  of  all  three  experiments  were  handed  to  you  in  my  office 
on  February  23,  1981. 

During  January,  1981,  the  surveying  firm  of  Harry  R. 

Feldman,  Inc.,  of  Boston,  Mass.,  used  conventional  surveying 
instruments  to  determine  the  same  baseline  vectors,  except  for  an 
overall  azimuth  rotation  which  remained  unknown.  After  I  had 

t 

given  you  our  radio  interferometry  results,  you  showed  me  for  the 
first  time  the  Feldman  survey  results.  We  compared  the  MITES 
with  the  Feldman  results  for  the  lengths  of  the  baselines,  but 
not  for  the  individual  horizontal  components  because  of  the  lack 
of  a  survey  determination  of  azimuth.  We  began  to  compare  the 
vertical  components  but  found  a  puzzling,  gr-oss,  discrepancy  that 
demanded  further  investigation. 

In  May,  1981,  I  found  the  apparent  explanation  of  the 
vertical  discrepancy:  one  simple  numerical  transcription  error 
by  Feldman.  In  July  1981,  we  obtained  the  result  of  an  accurate 
(±3.3  arcsec)  astronomical  determination  of  the  azimuth  of  one  of 
the  baseline  vectors  from  Bob  Magee  and  Mark  Nasson  of  CSDL 
(their  memo  no.  DSSG/M8 1 -1 9 ) .  I  have  now  combined  their  azimuth 
determination  with  Feldman's  original  between-basel ines  azimuth- 
difference  determinations  to  obtain  azimuths  for  all  baselines. 
Thus,  I  have  been  able  to  compute  all  three  rectangular 
components,  North  (N),  East  (E),  and  Up  (U) ,  of  each  baseline 
vector  according  to  the  combination  of  Feldman's  and  Maaee's 


results. 


In  the  present  memorandum  the  results  of  these  computations 
are  presented  and  compared  with  our  originally  reported 
interferometry  results  which  were  in  the  same  NEl)  coordinates. 

You  will  recall  that  in  order  to  tra  .sform  the  interferometry 
results  from  the  1950.0  mean  equator  al  coordinates,  in  which  our 
geometrical  calculations  had  been  pe  formed,  to  the  topocentric 
NEU  system,  I  assumed  that  the  U  axi  pointed  toward  oeocentric 
latitude  42°  37'  22",  longitude  71°  9'  19",  and  t.nat  the  E  axis 

was  parallel  to  the  equatorial  plane  of  date. 

Results 

In  Table  I  are  given  tie  N,  E,  I  components  and  the  length 
of  the  baseline  vector  that  extends  rom  our  laser-pad, 
manhole-rim,  X-cut,  temporary  survey  mark  (Feldman's  point  number 
4)  to  our  front  parking  lot,  screw-heac,  temporary  survey  mark 
(adjacent  to,  but  net  the  s  ime  as,  Feldman's  front  parking  lot 
"pin",  their  point  number  2'.  In  aciition  to  the  MITES  and  the 
survey  results  in  this  table  we  include  the  Geodimeter 
measurement  that  Feldman,  Inc.,  made  but  did  not  include  in  its 
final  report  to  you. 

In  Table  2  are  the  results  for  the  vector  from  ou  • 
screw-heac  temporary  surve  mark  in  the  oil-tank  area,  near  but 
not  the  same  as  Felcman's  lin  (their  point  no.  1),  to  our 
manhole-rim  X-cut  (Feldman  no.  4). 

In  Table  3  are  the  results  for  the  vector  from  our  frent-iot 


screw  head  to  our  oil-tank  screw  head. 
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In  Table  4  are  the  N  ,  E,  U  components  of  the  three  short 
offset  vectors  to  our  temporary  marks  from  the  nearby,  permanent 
marks:  one  from  the  NGS  disk  in  the  center  of  the  laser- pad 

manhole  (Feldman  pt.  no.  3)  to  the  rim  X-cut  (F.  pt.  no.  4);  one 
from  the  steel  pin  set  by  Feldman  in  the  front  lot  (F.  pt.  no.  2) 
to  our  nearby  screw  head;  and  one  from  the  steel  pin  set  by 

Feldman  in  the  oil  tank  area  (F.  pt .  no.  1)  to  our  nearby  screw 

« 

head.  The  N  and  E  component  values  given  in  Table  4  were 
computed  by  me  from  the  Feldman  and  Magee  distance  and  azimuth 
measurements.  The  U  component  of  the  offset  vector  was  not 
measured  by  the  surveyors  except  at  the  manhole;  I  measured  the 
two  others  myself  with  a  bubble  level  and  a  hand-held  scale. 

In  Table  5,  the  vector  misclosure  around  the  triangle  of 
baselines  is  shown  for  each  of  the  three  MITES  experiments.  (The 
conventional  survey  results  reported  by  Feldman  had  been  adjusted 
so  that  they  closed  exactly.  No  such  adjustment  was  made  for  the 
MITES  results. ) 

Discussion 

The  tables  indicate  that  both  the  MITES  and  the  conventional 
survey  results  are  accurate  at  about  the  estimated  level  of 
uncertainty,  about  5  mm  in  each  coordinate  except  for  the  second 
MITES  exD^riment  on  Dec.  29,  when  the  car  was  parked  next  to  one 
antenna.  The  uncertainty  for  this  experiment  is  higher  mainly 
because  of  the  poorer  satellite  geometry,  not  because  of  the  car. 
Note  that  the  baseline  of  Table  2  had  a  car  at  neither  end,  yet 
the  reoea tabi 1 i ty  of  its  vertical-component  determination  on 
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December  29  is  about  as  poor  as  for  the  two  baselines  which  had 
the  front-lot  antenna  in  common.  The  triangle  closure  error 
(Table  5)  is  extraordinary  large  for  the  Dec.  29  "car"  experiment 
because  in  this  experiment,  for  unknown  reasons,  each  of  the 
three  baselines  exhibited  a  centimeter-sized  upward  shift. 

Still,  it  seems  fair  to  say  that,  at  least  when  the  satellite 

geometry  is  good,  we  can  determine  all  three  components  of  a 

■ 

baseline  vector  that's  about  100  meters  long  within  about  5 
millimeters,  one  sigma,  by  MITES/GPS. 


Table  1.  North  (N),  East  ( E )■ ,  and  Up  (U)  components  and  length 
(L)  of  the  baseline  vector  extending  from  the  laser-pad  to  the 
front  parking  lot  temporary  survey  marks,  from  various 
determinations.  Uncertainties  are  subjectively  estimated,  and 
are  intended  to  represent  about  ±1  standard  deviation. 


Determination 


( Method/Date) 

MITES/Dec.  17 

MITES/Dec.  29 

MITES/Mean* 

MITES/Dec. 29/car** 

Survey 

Geod imeter 

Difference:  MITES/ 
Mean-Survey 

Difference:  MITES/ 
Me an -Geod imeter 


Baseline  Component 


U 


(mm) 

(mm) 

(mm) 

(mm) 

7 , 1 1 0± 

5 

-91,784+ 

5 

-1  ,289+ 

5 

92,068+  5 

7, 10  5± 

5 

-91 , 78  2± 

5 

-1 , 287± 

5 

9  2 , 065±  5 

7,107+ 

4 

-91,783+ 

4 

-1,288+ 

4 

92,067+  4 

7,107+ 

7 

-91,784+ 

7 

-1,278+ 

10 

92,067+  10 

7,1011 

5 

-91 ,780± 

5 

-1 ,291± 

5 

92,0631  5 

— 

— 

— 

92,074i 10 

6± 

6 

1 

u> 

1+ 

6 

3± 

6 

4  i  6 

— 

— 

— 

-7+11 

aver  age 

o  f 

the  Dec. 

17 

and  the 

first  Dec.  29 

result  shown,  although  computed  with  more  precision  before 
rounding  to  the  nearest  millimeter.  The  MITES/Dec . 29/ca r  result 
was  not  included  in  this  average. 

Obtained  from  observations  late  on  Dec.  29,  while  a  car  was 
parked  with  its  side  BO  cm  from  the  center  of  the  antenna  in  the 
front  parking  lot.  The  satellite  viewing  geometry  was 
s ion i  f l can tl v  poorer  during  these  observations  than  for  the 
earlier  ones;  thus,  the  baseline  uncertainty  is  magnified, 
especially  for  the  U  component. 


Table 


Z.  • 


North  (N),  East  (E).,  and  Up  (U)  components  and  length 


(L)  of  the  baseline  vector  extending  from  the  o  i  1- tank  to  the 
laser- pad  temporary  survey  marks,  from  various  determinations. 
Uncertainties  are  subjectively  estimated,  and  are  intended  to 
represent  about  -1  standard  deviation.  The  uncertainties  for  the 
survey  are  greater  for  this  baseline  than  for  the  others  because 
only  for  this  baseline  was  the  line  of  sight  obstructed  so  that  a 
dog-leg  had  to  be  surveyed. 


Determination  Baseline  Component 


N 

E 

U 

L 

(Method/Date) 

MITES/Dec.  17 

( mm ) 

-71 , 362± 

5 

( mm) 

101,1871 

5 

( mm) 

-5  61  5 

(mm) 

123,8201 

5 

MITES/Dec.  29 

-71 ,357± 

5 

101 ,1801 

5 

-521  5 

123,8111 

5 

MITES/Mean* 

-71,359± 

4 

101 ,184i 

4 

-54+  4 

123,816i 

4 

★  ★ 

MITES/Dec. 29/ca r 

-71 ,3651 

7 

101,1811 

7 

-44+10 

123,816110 

Survey 

-71,3471 

7 

101,1801 

7 

-4  7i  7 

123,8051 

7 

Difference:  MITES/ 

Mean-Survey 

-12± 

8 

4± 

8 

1 

1  + 

00 

111 

8 

★ 

Simple  arithmetic  average  of  the  Dec.  17  and  the  first  Dec.  29 
result  shown,  although  computed  with  more  precision  before 
rounding  tc  the  nearest  millimeter.  The  MITES/Dec . 29/ca r  result 
was  net  induced  in  this  average. 

★  *■ 

Obtained  from  observations  late  on  Dec.  29,  while  a  car  was 
parkec  with  its  side  60  cm  from  the  center  of  the  antenna  in  the 
front  parkinc  lot.  Note  that  neither  of  the  antennas  at  the  ends 
of  th<  baseline  vector  whose  components  are  aiver  in  this  table 
were  affected.  The  satellite  viewing  geometry  was  significantly 
poorer  curing  these  observations  than  for  the  earlier  ones;  thus, 
the  baseline  uncertainty  is  magnified,  especially  for  the  U 
component . 


Table  3.  North  (N),  East  (E)-,  and  Up  (U)  components  and  length 
(L)  of  the  baseline  vector  extending  from  the  front  parkina  lot 
to  the  oi I-tank  temporary  survey  marks,  from  various  determina¬ 
tions.  Uncertainties  are  subjectively  estimated,  and  are 
intended  to  represent  about  ±1  standard  deviation. 


Determination 

Basel ine 

Component 

N 

E 

U 

L 

(Method/Date ) 

( mm ) 

(mm) 

( mm ) 

(mm) 

MITES/Dec.  17 

64 ,242± 

5 

-9,4  00- 

5 

1  ,346± 

5 

64 ,940± 

5 

MITES/Dec.  29 

64 , 249± 

5 

-9,401+ 

5 

1  ,  342± 

5 

64,947+ 

5 

MITES/Mean* 

64,246+ 

4 

-9,400+ 

4 

•  1  ,344  + 

4 

64,944+ 

4 

MITES/Dec. 29/car** 

64,251  + 

7 

-9,406+ 

7 

1,365+  10 

64,950+ 

7 

Survey 

64,248+ 

4 

-9,399+ 

4 

1  ,33  8± 

4 

64,9^4+ 

a 

Difference;  MITES/ 
Mean-Survey 

-2± 

6 

-1± 

6 

6± 

6 

0± 

6 

* 

Simple  arithmetic  average  of  the  Dec.  17  and  the  first  Dec.  29 
result  shown ,  although  computed  with  mere  precision  before 
rounding  to  the  nearest  millimeter.  The  MI TES/Dec . 29/ca r  result 
was  not  included  in  this  average. 

*  -k 

Obtained  fren  observations  lat<=  on  Dec.  29,  while  a  car  was 
parked  with  its  side  60  cm  from  the  center  of  the  antenna  in  the 
front  parking  lot.  Th°  satellite  viewina  geometry  was 
significantly  poorer  durino  these  observations  than  for  the 
earlier  ones;  thus,  the  baseline  uncertainty  is  maanified, 
especially  for  tne  U  component. 


Table  4.  North  (N) ,  East  (E)  ,  and  Up  (U)  components  of  the  short 
offset  vectors  from  the  permanent  to  temporary  survey  marks. 
Estimated  uncertainty  is  ±1  mm  in  each  component  of  each  vector. 


Site 

Vector 

N 

E 

U 

From 

To 

( mm) 

(mm) 

(mm) 

La  se  r 

pad  manhole 

NGS  disk 

Rim  X- 

■cut 

262 

-5 

208 

Front 

lot 

Steel  pin 

Screw 

head 

-269 

27 

-16 

Oil  tank 

Steel  pin 

Screw 

head 

48 

227 

-1 

Table  5.  Vector  misclosure  around  the  triangle  of  baselines  for 
each  of  the  .^ITES  experiments.  The  entries  in  this  table  were 
obtained  by  summation  of  the  corresponding  entries  in  Tables  1, 

2,  and  3;  the  square  root  of  the  sum  of  the  squares  was  taken  for 
the  uncertainties. 


N 

E 

U 

Exper 

i  m  e  n  t 

( mm) 

(mm) 

(mm) 

Dec . 

27 

-10*  9 

4±  9 

1*  9 

Dec  . 

;  o 

-3*  9 

-3±  9 

3 1  o 

Dec  . 

+ 

19  /car 

-It  12 

-1 C±  12 

4  3±  17 

T 

Obtained  from-  observations  late  on  Dec.  29,  while  a  car  was 
pa  r  x ec  with  its  side  SO  cm  from  the  center  of  the  antenna  in  the 
front  oar  Kina  lot.  The  satellite  viewing  geometry  was 
significantly  poorer  during  these  observations  than  for  the 
earlier  ones;  thus,  the  baseline  uncertainty  is  maanified, 
especially  for  tne  U  component. 
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ACCURACIES  OF  BASELINE  DETERMINATIONS  BY  MITES 
ASSESSED  BY  COMPARISONS  WITH  TAPE,  THEODOLITE, 

AND  GEODIMETER  MEASUREMENTS 

C,  C.  Counselman  III,  S.  A.  Gourevitch,  R.  W. 
King,  T.  A.  Herring,  I.  I.  Shapiro  (Dept.  Earth 
&  Planetary  Sci.,  M.I.T.,  Cambridge,  MA  02139) 
R.L.  Greenspan  (C.S.  Draper  Lab.,  Camb.  ,MA  02139) 
A.E.E.  Rogers,  A.R.  Whitney,  R.J.  Cappalio  (NEROC 
Haystack  Observatory,  Westford,  MA  01886) 


On  Dec.  17,  1980,  portable  MITES  antennas  [ref. 
Bull.  Geod.  53,  139-163  (1979)]  were  set  atop 
three  survey  marks  near  the  Haystack  Observatory 
building.  These  antennas  yielded  radio  inter¬ 
ferometric  observations  of  the  NAVSTAR/GPS  satel¬ 
lites  which  were  analyzed  to  determine  the  vector 
baselines  between  the  survey  marks.  Observations 
were  repeated  12  days  later  with  different  anten¬ 
nas  on  the  marks.  For  the  determination  of  each 
baseline  on  each  day  all  of  the  observations  from 


the  entire  time  that  five  satellites  were  above 
20°  elevation — 1.3  hours — were  used.  No  data 
were  deleted  or  downweighted  and  the  same  para¬ 
meters  were  estimated  using  the  same  algorithm  in 
every  case.  On  each  day,  the  triangle  of  separ¬ 
ately  estimated  baselines  closed  within  1  cm  in 
each  vector  component;  and  for  each  baseline,  the 
two  determinations  agreed  within  1  cm.  In  Janu¬ 
ary,  1981,  the  firm  of  H.  Feldman,  Inc.,  surveyed 
the  triangle  conventionally  by  means  of  steel 
tape,  1"  theodolite,  laser  Geodimeter,  and  pre¬ 
cise  level.  The  MITES  experimenters  and  the 


conventional  surveyors  did  not  communicate,  but 
both  delivered  their  results  in  writing  to  a 
referee  (R.L.G.)  for  comparison.  The  results  for 
the  lengths  of  the  sides  of  the  triangle  were,  in 


millimeters : 

Tape 

MITES 

Geodimeter 


Side  1 
6^  ,9^  ±5 
6^.9^  -7 
n  .a . 


2 

92,063  ±5 
92,067  ±5 
92,07^10 


3 

123,805  ±5 
123,816  ±7 

n  .a . 


We  conclude  that,  at  least  for  short  baselines, 
the  resolution  of  interferometer  fringe  ambigui¬ 
ties  is  not  difficult;  also,  at  the  centimeter 


level  of  accuracy,  multipath  interference  is  not 
a  significant  problem  with  MITES. 


ACCURACY  OF  RELATIVE  POSITIONING  BY 
INTERFEROMETRY  WITH  RECONSTRUCTED 
CARRIER  GPS:  EXPERIMENTAL  RESULTS* 


Richard  L.  Greenspan 
Arthur  Y.  Na 
Joseph  M.  Prtvjemski 
James  D.  Veale 


The  Charles  Stark  Draoer  Laboratory,  Inc. 
Radio  Navigation  Division 
555  Technology  souare 
Cambridge,  MA  02139 


ABSTRACT 

Several  short  baselines  were  measured  by  means  of  interferometric 
(phase-difference)  processing  of  the  reconstructed  carrier  sianals  received 
from  four  GPS  satellites.  After  1  hour  of  observation,  the  baseline  esti¬ 
mates  agreed  with  surveyed  values  within  5  millimeters  in  each  coordinate. 
Ambiauities  in  the  carrier  phase  data  were  typically  resolved  within  the 
first  15  minutes  of  observation.  The  baseline  lengths  typically  converqed 
to  1  centimeter  of  their  final  value  within  1 5  to  30  minutes  of  observation. 
The  sensitivity  of  the  baseline  estimates  to  multipath  was  demonstrated  and 
theoretical  predictions  of  multioath  effects  were  confirmed.  This  experi¬ 
mental  prooram  demonstrated  the  feasibility  of  survevinc  short  baselines  to 
subcentimeter  accuracv  usino  highlv  portable  electronic  equipment.  Unpredic¬ 
table  components  of  propagation  delay  will  ->e  the  limiting  error  source  for 
longer  baselines. 


* 

This  work  was  sponsored  bv  the  U.s.  Air  porce  Geophysics  Laboratory  under 
Contract  Fi 9628-80-C0040. 
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INTRODUCTION 


The  Charles  Stark  Draper  Laboratory,  Inc.  (CSDL)  has  conducted  a  ser¬ 
ies  of  experiments  to  demonstrate  that  short  baselines  on  earth  can  be  meas¬ 
ured  in  maqnitude  and  direction  to  within  1  centimeter  or  less  of  error 
throuah  the  use  of  easilv  portable  electronic  equipment.  These  measurements 
were  obtained  by  interferometric  processinq  of  GPS  siqnals  received  at  two 
antennas  that  were  sited  at  the  ends  of  the  baselines  to  be  measured.*  The 
GPS  siqnals  were  amplified  and  cabled  to  a  GPS  receiver  wherein  the  narrow- 
band  carrier  siqnal  was  reconstructed  by  conventional  codetrack inq  proce¬ 
dures.  The  instantaneous  carrier  phase  for  each  satellite  in  view  was  sam¬ 
pled  periodically  to  form  the  set  of  observables  from  which  vector  baselines 
were  estimated.  The  followinq  sections  describe  the  experimental  setup,  sum¬ 
marize  the  data  collection  and  the  subsequent  results  of  baseline  estimation, 
and  qive  a  preliminary  interpretation  of  our  findinqs  with  respect  to  the  de¬ 
velopment  of  GPS  instrumentation  for  geodetic  applications. 

EXPERIMENTAL  DESIGN 


Fiaure  1  is  a  schematic  representation  of  the  hardware  configuration  of 
the  experiment.  The  antennas  (Miniature  Interferometer  Terminals  for  Earth 
Survevina  (MITES) )  and  preamplifiers  were  developed  for  CSDL  by  The  Stein- 
brecher  Corporation  in  collaboration  with  Professor  Charles  C.  Counselman  III 
of  MIT.  The  MITES  antennas  are  tuner  to  receive  the  GPS  L1  frequency 
1.57542  GHz.  The  antenna  is  functiorallv  comprised  of  two  orthogonal,  horizon¬ 
tal,  halfwave  resonant  dioole  elements  fed  in  Phase  quadrature  to  achieve 
circular  polarization  in  the  zenith  direction  (Reference  1).  The  dipoles 
are  mounted  on  a  cylindrical  brass  si  aft  at  a  height  of  3/8  wavelength  above 
a  metallic  horizontal  ground  plane,  1  sauare  yard  in  extent  (Figure  2).  The 
antenna  provides  a  nearlv  uniform  gatr.  from  the  zenith  to  about  a  20-degree 
elevation  above  the  horizontal.  Theiefore,  siqnals  from  all  satellites  at 
or  above  a  20-dearee  elevation  are  simultaneously  received  during  a  data 
collection  session  as  lonq  as  thev  ate  not  occulted  by  obstacles.  Signals 
from  individual  satellites  are  resol\  ed  in  the  GPS  receiver  by  means  of  their 
specific  code  structure. 

In  these  experiments,  the  antennas  are  precisely  located  over  survey 
markers  bv  means  of  a  pointed  brass  1  od  inserted  throuqh  a  vertical  hole 
bored  into  the  cylindrical  center  shift. 

The  antenna  output  is  coupled  to  the  preamplifier  unit  by  a  4-foot 
lenoth  of  RG-233  cable.  A  four-sect;on  cavity  bandpass  filter  with  25- 
meaahertz,  3-db  bandwidth  is  provider  to  reject  out-of-band  interference  with 
only  1  db  midband  insertion  loss.  A  limiter  provides  burn-out  protection  from 
strona  interfering  signals.  The  system  noise  fiaure  is  set  at  about  4  db 
bv  a  low-noise  GaAs  FET  preamplifier.  The  net  system  gain  is  approximately 


ror  the  purposes  of  this  paper,  shot t  baselines  are  those  for  which  prooa- 
oation  disturbances  are  highly  correlated  at  either  end  and  are  therefore 
cancelled  by  interferometric  proces:ina. 


M 


26  db.  The  preamplifier  output  is  connected  to  the  input  port  of  the  GPS 
receiver  by  approximately  35  meters  of  RG-214  double-shielded  coaxial  cable. 

The  GPS  receiver  is  the  breadboard  prototype  of  a  device  that  was  under 
development  at  CSDL.  The  hardware  consists  of  a  sinale  RF  channel  with  fre¬ 
quency  conversions  selected  to  receive  the  freauency.  For  this  experi¬ 
ment,  the  receiver  software  is  confioured  to  sequentially  track  four  satel¬ 
lite  sionals  from  each  baseline  antenna  after  they  have  been  acquired.  The 
receiver  implements  a  variation  of  the  Costas-loop  havinq  an  arctangent  error 
response  to  track  the  data-modulated  C/A  signals  received  from  a  single  sa¬ 
tellite  over  one  24-millisecond  interval.  At  the  end  of  each  interval,  it 
switches  under  software  control  to  process  the  next  satellite  in  the  next  in¬ 
terval.  Code  tracking  is  also  executed  to  make  GPS  timing  information  avail¬ 
able.  Data  recovery  to  obtain  satellite  ephemerides  is  performed  as  part  of 
the  initial  siqnal  acquisition  function. 

All  eiqht  receiver  channels  are  used  for  data  collection.  The  proce¬ 
dure  is  to  cable  antenna  outputs  from  both  ends  of  a  test  baseline  to  the  re¬ 

ceiver  and  to  process  up  to  four  satellites  from  each  end.  Thus,  the  cur¬ 
rent  configuration  implements  a  connected  element  interferometer.  However, 
we  encountered  one  source  of  Phase  drift  between  antenna  outputs  that  gives 
the  observables  the  characteristics  of  independent  element  interferometry. 

Sequencing  amonq  eight  channels  takes  approximately  192  milliseconds. 
Thus,  five  samples  of  the  carrier  phase  of  each  satellite  in  view  are  avail¬ 
able  every  second.  These  carrier  Phase  samples  are  ambiouous  by  multiples  of 

one-half  wavelenoth  because  of  the  arctangent  implementation.  However  the 
samplino  rate  is  fast  enouoh  that  carrier  cycles  are  not  lost.  Therefore,  a 
constant  time  offset  (equivalent  to  an  inteoral  multiple  of  half  wavelenaths) 
can  be  estimated  for  each  satellite  sianal  as  part  of  the  data  reduction. 


In  these  experiments  the  antenna/preamplif ier  units  were  located  on 
the  rooftop  at  CSDL  and  cabled  to  the  GPS  receiver  in  a  laboratory  one  floor 
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36-in.  SQUARE  REINFORCED  GROUND  PLANE 


a)  MECHANICAL  DESIGN 


b)  ELECTRICAL  DESIGN 

FIGURE  2.  MITES  ANTENNAS. 

below  the  roof  (Reference  2).  Fiqure  3  illustrate;;  the  orientation  of  eight 
markers  that  were  placed  on  the  rooftop  to  define  potential  baselines. 

These  markers  were  surveyed  by  taoe  and  by  theodolite  to  determine  the  base¬ 
line  lenaths,  the  elevation  of  the  markers  with  reSDect  to  marker  AT  and  the 
anales  of  the  baselines  with  respect  to  north.  The  anaular  measurements  are 


estimated  to  have  a  standard  deviation  of  3  arcseconds;  the  ooint  to  point 
distance  errors  are  estimated  to  have  a  sta  iriard  deviation  of  2  millimeters. 
Tables  1  and  2  summarize  the  survey  data. 


FIGURE  3.  SITING  OF  SURVEYED  EASELINES  ON  CSDL  ROOF. 

During  anv  data  collection  session,  two  antennas  are  centered  over  the 
surveyed  points;  thev  are  leveled  and  oriented  within  a  few  degrees  of  north 
using  a  bubble  level  and  compass.  The  orientation  is  important  because  the 
antenna  patterns  will  introduce  an  electrical  phase  angle  epual  to  their  dif¬ 
ference  in  azimuth  orientations.  Because  the  two  antennas  are  essentially 
identical  with  a  circular  svmmetry,  their  phase  centers  can  be  taken  at  the 
height  of  the  around  plane  directly  above  the  survey  markers.  1t>e  antennas 
have  screw-tvoe  levelino  leas;  therefore  one  must  measure  the  height  of  each 
antenna  above  the  marker  in  order  to  properly  compare  the  interferometry  re¬ 
sults  with  the  conventional  survey.  This  was  done  by  measuring  the  heiaht  of 
the  central  shaft  above  the  marker  plate  for  each  experiment.  This  measure¬ 
ment  is  accurate  to  about  1  millimeter. 


All  experiments  were  conducted  when  three  or  more  GPS  satellites  were 
visible  above  a  20-dearee  elevation  as  seen  from  Cambridge,  Massachusetts. 


This  limited  the  potential  viewinq  time  to  about  78  minutes  each  day  of  which 
about  16  minutes  were  required  for  siqnal  acquisition,  eohemeris  collection, 
and  determination  of  GPS  svstem  time.  Therefore,  data  collection  typically 
scanned  about  1  hour.  As  a  result,  only  one  baseline  was  measured  per  day. 
Fiqure  4  (Courtesy  of  C.  Counselman)  illustrates  the  azimuth  and  elevation 
histories  of  the  visible  satellites  durinq  a  data  collection  session.  Since 
the  GPS  satellites  are  in  a  nearly  12-hour  orbit,  the  same  pattern  repeats 
from  day  to  day  at  or oar essively  earlier  times  (4  minutes  per  dav) . 

TABLE  1 .  HORIZONTAL  COMPONENTS  OF  REFERENCE  BASELINES . 


DISTANCE 

(meters) 

ANGLE  FROM  NORTH 
(degrees) 

A1A4 

19.6470 

44.7808 

A4A2 

40.9507 

110.9572 

A1A2 

52.0859 

90.7712 

A2A7 

31.6031 

148.5794 

A3A7 

20.9295 

48.4691 

A8A3 

22.0182 

131.4213 

A5A8 

29.8460 

109.2713 

A5A3 

50.9205 

118.6541 

A1A5 

19.9920 

154.3862 

A5A6 

8.2000 

22.721 1 

A1A6 

14.8620 

130.0446 

A5A4 

31 .5772 

10.2689 

A6A4 

23.6363 

5.9788 

TABLE  2.  ELEVATION  OF  REFERENCE  MARKERS. 


ELEVATION 

(cm) 

A1 

0.00 

A2 

-11.45 

A3 

-10.65 

A4 

-8.80 

A5 

-8.9  5 

A6 

-8.10 

A7 

-0.60 

A8 

-0.3  5 

EXPERIMENTAL  OBSERVATION 

Data  collection  beaan  on  Auoust  13,  1981  with  observations  of  baseline 
A 1 A2 .  The  first  surveys  used  carrier  phase  differences  (inteqrated  Doppler) 
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-  CIRCLED  NUMBERS  ARE  GPS  NAVSTAR  SATELLITE  DESIGNATORS 

-  SATELLITE  TRACKS  ARE  DARKENED  TO  INDICATE  SIMULTANEOUS 
VISIBILITY  ABOVE  20*  ELEVATION 

FIGURE  4.  VISIBILITY  OF  GPS  SATELLITES. 

as  the  observable.  The  resultinc  baselines  disaoreed  from  the  survey  bv  sev¬ 
eral  centimeters.  Althouch  this  was  encouraainc,  we  beaan  an  intensive  re¬ 
view  of  the  data  collection  orocedures  and  the  data  reduction  software  in  or¬ 
der  to  determine  the  cause  of  unanticipated  trends  in  the  oost-fit  residuals. 
These  trends  were  the  evidence  that  led  to  the  discoverv  of  a  sianificant 
timinc  error  in  the  processing  software.  By  the  end  of  SeDtember  1981,  the 
orocessina  software  was  debuaaed  and  extended  to  include  the  option  for  phase 


difference  processinq,  as  well  as  integrated  Doppler  difference  processing. 
Data  collection  for  baseline  measur ements  was  concentrated  in  October  1981, 
with  further  observation  sessions  durinq  December  to  assist  in  diagnosing  the 
measurement  error  sources.  Table  3  summarizes  the  data  base  collected  during 
this  period.  Weather  conditions  durinq  baseline  measurements  ranged  from 
clear,  calm,  and  dry  with  tenroeratures  around  22  C  to  wet  and  windy  with  tem¬ 
peratures  around  10  C. 


TABLE  3.  BASELINE  MEASUREMENT  DATA  3ASE. 


Dates 

Baselines 

Comments 

October 

1,3 

Null 

Used  to  estimate  noise  level  attributable 
to  receiver  electronics 

October 

7,8,9,13 

A5A4  \ 

Used  to  estimate  measurement  repeatability 

October 

14,15 

*,AS  , 

and  closure  of  trianqular  survey 

October 

16 

A1A4  ' 

October 

21,23,26 

A5A4 

Diagnostic  tests  for  observed  phase  drifts 

October 

27 

A5A4 

Reflecting  mirror  placed  to  produce  multipath 

December 

15,16,17 

Null 

Diagnostic  for  contribution  of  electronics  to 
oost-fit  residuals 

DATA  REDUCTION 


The  GPP  receiver  generates  sets  of  eight  phase  measurements  taken  se¬ 
quentially  over  a  192-millisecond  interval  which  corresponds  to  the  samplinq 
oer iod  of  the  receiver's  ohase  tracking  filters.  For  a  baseline  AB,  the  ob¬ 
servations  are  sequenced  as  Satellite  n  at  A,  Satellite  n  at  B,  Satellite 
n  +  1  at  A,  and  so  forth.  The  receiver  software  also  outputs  estimates  of 
the  carrier  ohase  rates  and  accelerations  so  that  all  observables  can  be 
referred  to  a  common  epoch.  Experience  has  shown  that  the  ohase  rate  terms 
are  adeauate  for  this  interpolation;  there  is  no  aoparent  improvement  of  the 
post-fit  residuals  when  the  acceleration  term  is  included. 

Th«  bandwidth  of  the  carrier  phase-tracking  filter  is  aDproximately 
1/2  hertz.  Therefore,  the  minimum  period  between  useful  sets  of  observables 
ir  approximately  2  seconds.  Although  all  of  this  data  are  available,  we  have 
t.vpicaliv  chosen  to  process  data  at  19.2-second  intervals  as  a  qood  compro¬ 
mise  between  baseline  measurement  accuracy  and  processing  burden.  We  will 
discuss  the  dependence  of  measurement  accuracy  on  the  number  of  observations 
.  r>  the  following  section. 

The  basic  set  of  observables  is  the  difference  in  carrier  phase  for 
recension  of  each  satellite  at  each  end  of  the  baseline.  These  differ- 
are  formed  after  the  sequential  observations  are  interpolated  to  a  com- 
-  •  epoch.  These  observables  were  processed  in  three  steps.  In  the  first 
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step,  a  linear  least  squares  algorithm  with  equal  weighting  of  residuals  was 
applied  to  generate  a  trial  solution  for  the  ambiguous  baseline  components. 
The  following  set  of  equations  was  fit  to  the  data  to  estimate  eight  param¬ 
eters,  including  three  baseline  components,  frequency  offset,  and  one  phase 
offset  for  each  satellite. 


j  =  b  ‘  sij  +  a  +  i (66) 


j  =  1,2, 3, 4 


where 


t.  .  =  the  interferometric  phase  difference  for  the  i  observa- 

1')  tion  of  satellite  j 

■  ^ 

b  =  the  unknown  baseline  vector 

-► 

=  a  unit  vector  in  the  direction  of  satellite  j  at  time  i 

=  an  unknown  phase  offset  for  observations  of  the  j*”*1 
satellite 

6e  =  the  phase  shift  produced  by  an  unknown  frequency  offset  be¬ 
tween  the  receiving  systems 

The  resulting  values  of  a  were  converted  to  integer  multiples  of  A/2,  and  a 

limited  search  around  these  integers  was  executed  to  determine  the  set  of  inte 
gers  that  minimizes  an  error  criterion.  Call  these  integers  N  .  The  final 

step  of  the  processing  is  to  refit  the  baseline  by  a  five-parameter  least 
squares  procedure  giving  the  solution  to  the  equations 

~  A/2  =  b  •  s^  +  i(60)  +6  3  1,2, 3, 4 

where  6  is  included  to  fit  cable  length  differences  that  are  a  fraction  of 


Table  4  summarizes  the  results  of  the  baseline  estimation.  The  base¬ 
line  components  are  presented  in  local-level  (north-east-up)  coordinates. 
There  was  no  data  editing  in  these  runs,  although  subsequent  analysis  sug¬ 
gests  that  a  short  interval  of  data  taken  immediately  after  initial  satellite 
acquisition  should  be  suppressed  to  avoid  contamination  by  the  transient  re¬ 
sponse  of  the  carrier  tracking  filters.  Each  baseline  listed  in  Table  4  was 
measured  on  a  different  day;  one  antenna  was  moved,  repositioned  and  leveled 
from  day  to  day.  The  largest  disagreement  between  the  conventional  survey 
and  the  phase  processed  data  is  5.7  millimeters  (east  component  of  A5A4) ,  and 
five  of  the  nine  differences  are  less  than  0.4  millimeters.  The  rms  post- 
fit  residual  for  phase  estimation  is  about  2  millimeters.  The  horizontal 
disagreements  for  Doppler  processing  range  between  1  to  4  centimeters,  with 
larger  errors  in  the  vertical  component,  reaching  12  centimeters  in  one  case. 
Vector  misclosure  around  the  triangle  of  baselines  is  -2.2,  -2.2,  -0.3  (mm) 
in  the  coordinate  directions  (N ,  E,  U) ,  respectively. 


TABLE  4.  EASELIN’E  SURVEYS. 


Basel  me 

Survey 

Mode 

North 

tmeter s ) 

East 

(meters) 

up13' 

(meters) 

RMS  Residual t4) 
(millimeters) 

A  1  A4 

Convent lonal 

13.9453 

13.8390 

-0.0880 

— 

A1A4 

Inter  ferometric 

(10/16/81) 

i 1 ) 

Phase  , 

_  .  (2) 
Doppler 

13.9481 

13.9581 

13. 8388 
13.8259 

-0.0883 

-0.1256 

1.99 

1.39 

A 1  AS 

Conventional 

-17.1253 

8.2101 

-0.0895 

— 

A  1 A5 

Inter  ferometric 

(10/15/81) 

Phase ( 15 
Doppler 

-17.1256 

-17.1477 

8.2131 

8.2471 

-0.0878 

-0.0111 

2.07 

1.47 

A5A4 

Conventional 

31.0714 

5.6292 

0.0015 

— 

A5A4 

Interferometric 

(10/13/81) 

I 

Doppler 

31.0715 

31 . 0981 

5.6235 

5. 5923 

-0.0008 

-0.1217 

2.02 

1.50 

- 

NOTES: 

i 

|  (1)  Sinqle-d  i f f erencr-d  pnase  processing. 

j  (2)  Sinqle-d if f er enced  integrated  Doppler  processing. 

|  (3)  Conventional  values  of 
feromevnc  measurements 
i  above  the  survey  marxs. 

the  ’Up"  component  are  from  Table 
are  adjusted  for  the  heights  of 

2;  the  inter- 
the  antennas 

t 

(  (4)  187  data  sets  spaced  at 
i  vat  ions  are  fit  to  five 
;  parameters. 

1  R  .  2  2 -second  intervals  over  1  hour.  Phase  obser- 
parameters  and  Doppler  observations  fit  to  three 

Table  5  summarizes  the  results  of  repeated  measurements  of  baseline 
A5A4.  The  baseline  components  typically  differ  by  less  than  2  millimeters, 
except  for  one  measurement  of  the  vertical  component  which  differs  by  about 
(  millimeters.  The  standard  deviations  are  0.76,  1.60,  2.40  (mm)  in  the 
coordinate  directions  (N,  F  U) ,  respectively. 

DISCUS5ION 

The  results  presented  in  Table  4  are  repeatable.  This  experiment 
gives  primary  evidence  for  the  following  conclusions: 

(It  Interferometric  processina  of  carrier  phase  differences  ob¬ 
served  in  reception  of  GPS  signals  can  be  used  to  measure 
short  baselines  to  an  accuracy  of  better  than  1  centimeter  in 
length  and  in  each  coordinate. 


TABLE  5.  REPEATABILITY  OF  BASELINE  ESTIMATES 
(Baseline  A5A4) 


Date 

North  (mi 

East  (m) 

Ud  (m) 

l 

October  7 

-31.0731 

-5.6226 

0.0253 

8 

-31.0735 

-5.6234 

0.0306 

9 

-31.0730 

-5.6227 

0.0312 

_ 1L_J 

-31.0715 

-5.6235 

0.0313 

(2)  For  lonqer  baselines,  the  primary  error  source  will  be  sta¬ 
tistically  independent  propagation  delays  along  the  path  from 
the  satellite  to  each  baseline  antenna.  These  propaaation  re¬ 
siduals  are  now  believed  to  contribute  about  1  to  2  centi¬ 
meters  error  to  lonq  baseline  measurements  (References  3 

and  4).  For  baselines  of  1000  kilometers  or  more,  satellite 
ephemeris  errors  could  become  siqnificant. 

(3)  The  quoted  accuracy  from  GPS  observations  can  be  achieved  us¬ 
ing  hiqhly  portable  equipment  includinq  a  ruaqed,  low-qain  an¬ 
tenna  (e.q.,  MITES),  a  sincle-channel  rapidly  sequencing  GPS 
receiver  and  data  acouisition  system  under  control  of  a  small 
minicomputer,  and  a  high  qualitv  crystal  oscillator  havinq  a 

short-term  stability  on  the  order  of  1  part  in  1011  over  100 
seconds . 

We  executed  several  variations  on  the  basic  data  collection  and  re¬ 
duction  procedures  in  order  to  sharoen  these  conclusions.  The  first  varia¬ 
tion  involves  measurements  of  a  "null"  baseline.  This  is  a  procedure  where¬ 
in  the  signals  from  one  antenna  are  split  into  two  paths  and  input  to  the 
GPS  receiver  as  if  they  came  from  two  different  antennas.  It  is  used  to  es¬ 
timate  the  noise  level  of  the  baseline  measurement  procedure. 

Figure  5  illustrates  an  overlay  of  the  Phase  residuals  for  all  satel¬ 
lites  observed  durino  the  null-baseline  test  of  October  1,  1981.  In  this 
test,  the  signal  splitting  occurred  at  the  input  to  the  G*>S  receiver.  Thus 
the  antenna,  preamplifier,  and  all  cabling  is  common  to  both  ends  of  the 
"baseline".  Data  were  collected  for  35  minutes;  for  this  data  set,  the  least 
squares  fit  for  the  baseline  components  yields  1.8,  -0.3,  -0.8  millimeters 
with  an  rms  oost-fit  residual  of  0.96  millimeters.  The  results  converge  to 
this  level  within  15  minutes  of  observation.  The  residuals  appear  to  be  ran- 
omlv  distributed  with  no  discernable  trends.  We  conclude  that  the  sequencing 
amonc  receiver  channels  is  fast  enough  that  receiver  electronics,  includinq 
phase  instability  of  the  receiver  frequency  standard  (FTS-Model  1000),  is  a 
necliqible  error  source. 

The  rms  oost-fit  residual  converges  to  the  standard  deviation  of  the 
observation  noise  when  that  noise  is  white.  The  rms  residual  can  be  predicted 
from  SNP  considerations  as  follows.  The  car r ler-to-noise  ratio  input  to 
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FIGURE  5.  POST-FIT  PHASE  RESIDUALS  FOR  "ZERO- BASELINE" 
OBSERVATION  OF  OCTOBER  1,  1981. 


the  sequential  trackinq  filters  was  estimated  at  35  +  2  db-Hz  and  the  track- 
inq  looo  noise  bandwidth  was  approximately  1/2  Hz.  For  C/n  =  35  dB-Hz, 
the  rms  error  in  measurinq  phase  distances  is  given  by 


oX 


/2  • 


B  (Hz) 
n 


2tt  C/n  (dB-Hz) 
o 


0. 55  mm 


which  is  within  a  factor  of  two  of  the  observed  phase  error. 

""he  zero-baseline  data  was  also  teorocessed  in  the  "double-difference’ 

mode 


Observable  =  (^.  -  -  (*2j  -  *2k  ) 

where  .  .  is  the  carrier  chase  of  satellite  j  as  viewed  viewed  at  the  it^1 
antenna.  The  rms  residual  increases  bv  almost  -J2  to  1.33,  which  is  consis¬ 
tent  with  the  assumption  that  these  residuals  are  uncorrelated. 

Fiqure  6  illustrates  the  phase  residuals  for  the  first  processinq  of 
the  observations  of  October  13,  15,  16  cited  in  Table  4.  There  is  clear  evi¬ 
dence  of  a  chase  slooe  in  the  data.  The  scatter  around  this  slope  appears  to 
have  a  neak-to-peak  ranqe  about  3  millimeters.  The  siqnificant  feature  of 
Fiqure  6  is  that  the  trend  is  common  to  all  satellites.  That  rules  out 
satellite-dependent  antenna  or  multipath  effects  and  suqaests  that  there  is  a 
small  time-varvinq  phase  drift  between  the  two  antenna  channels.  The  obser¬ 
vation  of  this  slow  drift  term  was  the  motivation  to  include  a  phase  drift 
(frequency  offset)  parameter  in  the  baseline  estimation  algorithm.  The  re¬ 
sults  of  Table  4  were  produced  with  the  Phase  slope  estimated.  Fiqure  7  il¬ 
lustrates  the  post-fit  residuals  usina  this  improved  estimator. 

Further  tests  were  conducted  in  December  to  diaanose  the  source  of 
the  drift.  "Tie  results  of  the  first  null-baseline  test  exclude  the  receiver 
as  the  source  of  the  drift.  The  remaininc  candidates  were  the  preamplifier 
units  or,  perhaps,  the  cavity  filters  or  cables. 
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(b)  BASELINE  15,  OCTOBER  15,  1981 
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(c)  BASELINE  14,  OCTOBER  16,  1981 


FIGURE  6.  SMOOTHED  PC3T-FIT  PHASE  RESIDUALS  FOR  MEASURED 


(b)  BASELINE  15,  OCTOBER  15,  1981 


FIGURE  7.  SMOOTHED  POST-FIT  PHASE  RESIDUALS  FOR  MEASURED 
BASELINES:  PHASE  DIFFERENCES  SLOPE  ESTIMATED. 
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In  the  zero-baseline  test  of  December  IS,  the  antenna  was  cabled  to 
one  preamolif ier  whose  outDut  was  sol  it  and  cabled  to  the  receiver.  This  is 
equivalent  to  the  first  zero-length  baseline  experiment,  except  that  there 
is  an  opoortunity  for  the  cables  to  be  a  source  of  phase  error.  Fiqure  8 
illustrates  the  phase  residuals  on  a  sample  by  sample  basis  and  with  smooth¬ 
ing  over  5  adjacent  samples  (96  seconds) .  There  appears  to  be  a  slope  of 
about  4  millimeters  per  hour.  Since  the  preamplifiers  are  common  to  both 
paths,  this  slope  suqqests  that  there  was  a  change  in  the  length  of  the  ca¬ 
bles.  The  thermal  coefficient  for  electrical  length  of  a  polyethylene  insu¬ 
lated  cable  (RG-214)  ranoes  from  100  to  250  ppm /‘C  (Reference  4).  For  a  tem¬ 
perature  rise  of  2‘C/hour  (during  mid-day  December) ,  we  could  expect  that 
phase  slopes  on  the  order  of  5  millimeters  per  hour  or  more  could  be  attrib¬ 
uted  to  the  20-meter  difference  in  the  lenaths  of  cables  that  were  exposed  to 
the  chanainq  thermal  environment.  The  experiment  was  repeated  with  the  an¬ 
tenna  output  split  and  amplified  separately  before  beinq  cabled  to  the  re¬ 
ceiver.  A  phase  slope  of  2  or  4  millimeters  per  hour  was  observed  on  suc¬ 
cessive  davs  when  different  pairs  of  amplifiers  were  used.  The  use  of  inde¬ 
pendent  amplifiers  did  not  increase  the  phase  slope.  This  confirms  our  sus¬ 
picion  that  differential  expansion  of  the  cables  produced  the  observed  ef¬ 
fect. 
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Another  question  we  considered  was  the  rate  at  which  the  baseline  es¬ 
timates  converqe,  and  in  particular,  how  much  time  is  reouired  to  resolve 
the  ambiauities  in  carrier  ohase  differences.  The  result  of  more  than  ten 
baseline  determinations  is  that,  in  most  cases,  the  ambiguities  were  correct - 
lv  resolved  for  the  trial  baselines  from  10  to  15  minutes  of  observation,  as 
lonq  as  there  were  several  data  Doints  per  minute.  The  exceptions  occurred 
when  the  phase  slope  (that  we  attribute  to  cable  stretch)  was  large.  In 
these  cases,  as  much  as  35  minutes  was  required  to  properly  resolve  ambiqui- 
ties.  In  most  other  cases,  the  baseline  lenqth  was  estimated  with  less  than 
1  centimeter  error  within  10  to  15  minutes.  After  that,  the  errors  steadilv 
decreased  to  the  level  of  a  few  millimeters.  We  conclude  that: 

(1)  The  rms  residual  converqes  rapidly  after  ten  or  more  data 
points  are  processed. 

(2)  The  accuracy  of  the  baseline  estimate  improves  with  the  span 
of  the  observations  and  the  number  of  data  points.  If  the 
number  of  data  points  is  fixed,  then  errors  decrease  if  the  ob¬ 
servation  period  is  increased. 

The  final  question  we  address  is  the  influence  of  multipath  on  the 
baseline  estimation.  Baseline  A5A4  w?s  remeasured  on  October  26  and  then  on 
October  27  a  4-foot  hiqh,  2-foot  wide,  1/8  inch  aluminum  plate  facing  northwest 
was  mounted  vertically  behind  the  antenna  at  A5,  as  illustrated  by  Figure  9. 

The  procedure  outlined  in  Reference  5  was  used  to  predict  the  incidence  of 
multioath  during  a  1-hour  observation  session.  The  theory  predicts  that 
siqnals  from  satellites  labeled  5,  6,  and  8  would  not  be  reflected  in  the 
mirror;  siqnals  from  satellite  4  'would  reflect  for  a  short  interval  but  these 
reflections  would  not  produce  significant  phase  deviations;  siqnals  from 
satellite  R  would  be  reflected  durinq  roost  of  the  observation  period  and 
this  multiDath  would  cause  a  peak-to-oeak  phase  offset  of  about  19  milli¬ 
meters.  Figure  10  illurtrates  the  observed  phase  residuals .The  upper  plot 
shows  the  residuals  for  October  26  (no  plate)  and  residuals  for  October  27. 

In  both  cases,  the  residuals  have  been  smoothed  over  five  successive  samples 
(19.22  seconds  per  sample)  to  suppress  short-term  noise  fluctuations.  Com¬ 
parison  of  these  residuals  shows  that  the  theory  is  qualitatively  accurate 
and  remarkably  close  to  beinq  auantit  itive.ly  accurate.  The  N,  E,  U  compo¬ 
nents  of  baseline  A5A4  estimated  from  the  data  of  October  26  are  -31.0423, 
-5.6371,  0.0003  (meters),  resoectivel".  The  baseline  length  is  31.5500 
meters.  These  results  deviate  from  prior  measurements  of  A5A4  because  the 
antenna  was  not  accurately  sited  over  the  A5  marker  on  these  two  days. 

The  observations  of  October  27  yield  baseline  components  of  -31.0048,  -5.6691, 
and  0.1620  (meters).  The  baseline  lenqth  is  estimated  as  31.5792  meters. 

This  inauspicious  placement  of  the  reflecting  plane  demonstrates  the  possi¬ 
bility  that  multioath  reflectina  into  a  poorly  sited  antenna  could  produce 
baseline  measurement  error  of  a  few  centimeters. 

CONCLUSIONS 

We  have  demonstrated  that  reconstructed  carrier  processing  of  GPS  sig¬ 
nals  can  be  used  to  achieve  subcentimeter  accuracy  in  the  measurement  of 
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(a)  PHASE  RESIDUALS  FOR  BASELINE  54,  OCTOBER  26.  1981 


(b)  PHASE  RESIDUAL  FOR  BASELINE  54  WITH  REFLECTING  PLATE,  OCTOBER  27.  1981 

SATELLITE  IDENTIFIER  4  WWW 
6  ////// 

8 - 

9  min 


FIGURE  10.  INFLUENCE  OF  MULTIPATH  ON  POST-FIT  PHASE 
RESIDUALS. 

short  baselines.  With  the  advent  of  compact  GPS  receivers,  the  field  in¬ 
strumentation  deluding  antennas,  clocks,  receivers,  and  data  collection  will 
be  hiqhlv  portable.  The  reauired  GPS  technology  involves  carrier  trackina 
and  rapid  sequencing  of  C/A  code-tracking  Ioods.  Second— gener at  ion  portable 
GPS  equipment  embodying  these  functions  is  currently  under  development  by 
the  Defense  Deoartment. 


We  believe  that  the  Drimarv  error  source  for  longer  baselines  will  be 
uncor rectable  propagation  delays  and  satellite  ephemeris  errors.  Multipath, 
cable  stretch,  clock  drift,  time  synchronization,  and  thermal  noise  will  not 


be  sianificant  error  sources  in  a  well-desiqned  and  well-sited  system;  nor 
will  the  suppression  of  these  effects  incur  any  sianificant  penalty  in  cost 
or  convenience  to  Dotential  users. 
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Introduction 

Multipath  is  a  potential  source  of  errors  in  measuring  baselines  by 
means  of  carrier  phase  interferometry.  Reflections  of  the  radio  signal 
off  of  nearby  objects  can  combine  with  the  direct  path  signal  to  distort 
the  measured  carrier  pnase.  This  multipath  distortion  is  not,  in 
general,  common  to  the  two  ends  of  the  baseline  because  the  topography 
differs,  and  is  not  common  to  two  satellites  tracked  simultaneously  at 
one  end  because  the  satellites  are  at  different  points  in  the  sky. 
Formation  of  the  interferometric  phase  observable  with  a  single 
difference  (phase  at  one  end  subtracted  from  phase  at  the  other  end  for 
simultaneous  measurements)  cancels  errors  that  are  common  to 
measurements  made  at  the  two  ends  of  the  baseline.  Use  of  a  second 
difference  (phase  due  to  one  satellite  subtracted  from  phase  due  to 
another  satellite)  cancels  errors  that  are  common  to  measurements  made 
at  either  end  which  are  the  same  for  more  than  one  satellite  —  examples 
of  such  errors  include  station  clock  errors  and  cable  stretch  due  to 
thermal  loading.  Multipath,  however,  falls  into  neither  of  these  two 
classes  of  "removable”  errors.  It  is  therefore  of  interest  to  learn  how 
large  the  multipath  distortion  may  be  and  in  what  situations  may 
multipath  be  a  significant  error  source. 

We  derive  an  analytic  expression  for  the  contribution  of  multipath  to 
the  received  pnase  for  the  specific  case  of  a  vertical  plane  reflector. 
The  intent  is  to  give  ourselves  a  qualitative,  and  to  some  extent, 
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quantitative  feel  for  the  size  of  the  multipath  distortion.  This  result 
is  a  special  case  of  a  general  expression  which  we  give  in  operator 
notation.  The  general  expression  covers  antennas  like  the  HITES  [lj 
antenna,  chat  is,  a  quadrature  fed  pair  of  dipoles  mounted  at  right 
angles  to  each  other  in  che  horizontal  plane  a  distance  g  (for  HITES,  g 
■  3  X / 8)  above  a  horizontal  ground  plane. 

An  experiment  has  been  conducted  on  the  CSDL  roof  which  demonstrates  the 
observability  of  multipath  due  to  a  rectangular  reflector  mounted  is  the 
vertical  plane  a  short  distance  from  one  antenna  of  a  baseline.  The 
residuals  after  the  fit  resemble  the  theoretical  phase  error  fairly 
closely.  The  experiment  is  described  and  qualitative  results  are 
presented. 

Theoretical  Development 

Tne  following  analysis  uses  the  operator  notation  and  analysis 
techniques  employed  by  Counselman  and  Gourevitch  [2J.  The  situation  is 
depicted  in  Figure  1.  Consider  a  baseline  consisting  of  antennas  A  and 
B,  each  with  its  own  receiver.  Only  antenna  A  has  reflecting  objects 
around  it;  antenna  B  is  assumed  to  be  on  level  ground,  which  is  by 
assumption  an  extension  of  the  antenna  ground  plane.  A  satellite  at 
zenith  angle  6  and  east  azimuth  $  (referenced  to  north)  Is  observed  by 
both  receivers,  which  track  the  phase  of  the  received  carrier.  The 
observable  is  the  difference  in  received  phase  sampled  simultaneously, 
or  nearly  so.  It  has  been  shown  (3)  that  in  the  absence  of  errors,  the 
phase-connected,  ambiguous  observable  is  equal  to  ♦ftji  where  r. 

denotes  range  to  the  satellite  from  receiver  i,  and  the  term  nit 
accounts  for  half  wavelength  ambiguities  which  must  be  removed. 

The  reflection  produces  a  relatively  delayed  signal  that  interferes 
coherently  with  the  desired  direct-path  signal.  The  reflected  signal 
arrives  attenuated  in  strength  and  shifted  in  phase  with  respect  to  the 
direct  signal  from  the  satellite  because  of  the  extra  propagation  delay, 
azimuth  dependence  of  the  antenna  pattern  phase  characteristic,  and  the 
reversal  of  sense  of  circular  polarization  due  to  reflection.  The 
effect  on  the  measured  signal  can  be  seen  qualitatively  in  the  phasor 
diagram  in  Figure  2.  In  the  diagram  c  s  represents  the  portion  of  the 
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signal  received  : roc  the  satellite  anc  :  r  represents  tne  portion  of 
the  signal  received  f roc  the  reflector.  A>  tne  satellite-anteana- 
rtflector  geoaetry  changes  with  tiae.  the  enepoiat  c:  ;r  traces  out  a 
locus  that  is  not  necessarily  circular.  Indeed,  given  the  maerout 
factors  contributing  tc  received  aultipath  phase,  the  length  end 
orientation  of  c •  with  respect  to  e»  eight  well  look  arbitrary.  Uc 
wish  tc  itnos.*  nor*  oetails  concerning  the  relationship  between  4  r  and 
}tl  reasonable  limits  on  the  phase  error  tt ,  and  rate  of  variation 
of  «(  as  the  satellite  moves.  One  of  the  principal  goals  of  this 
section  is  to  characterize  received  aultipath  phar  error  by  deriving  an 
expression  for  the  special  caee  of  a  vertical  reflector. 

In  the  notation  of  [  the  direct  satellite  signal  received  by  the 
quadrature  fed  crossed  dipole  antenna  is  given  by  the  following 
expression,  neglecting  an  asplitude  factor  chat  has  been  set  to  unity 
throughout  this  discussion: 

A*  <?•!*>■'!  p>  »<*U*l3>)  "> 

where  j  ’•  ^  is  the  complex  colmn  vector  representing  the  orientation  of 
the  quadrature-fed  horizontally  mounted  crossed  dipoles  of  the  KITES 
antenna: 


v  VT  • 


^  pw]  is  the  complex  conjugate  transpose  of  [p}  ;  D  is  an  identity 
operator  that  performs  a  rotation  in  azimuth  and  zenith  angle  to  locate 


the  transmitter  in  space: 


f  t#t0 


cc*4> 


*i'i»0tel$ 


r\  m 


I  W (  u  the  propagation  (row)  vector  of  the  si gnai  it oc  its  source  and 
nae  magnitude  Sf/i;|g^r*  tho  ground  plan*  vector  auditors  the  range 
vector  free  me  antenna  to  the  satellite.  The  first  element  ir.  the** 
vectors  is  the  *  or  north  component ,  th«  second  is  the  y  or  west 
component ,  and  the  third  is  the  s  or  up  component.  '-her.  the  quantities 
listed  abov*  at#  substituted  into  (1 ‘  w#  get  t.»*  following  useiul 
expression: 

A  *  |(l  tevfi)  € 

The  bracketed  tens  ir.  (?)  is  the  antenns  pattern  for  the  received 
electric  field.  The  expression  m  (?)  gives  the  complex  amplitude  of  a 
received  signal  due  tc  a  source  a:  rang*  r,  tentth  angle  f ,  aciauth  s, 
transcutting  circularly  polarized  radiation  in  the  direction  of  the 
artenna. 


A  similar  expression  can  be  derived  for  signals  reflected  into  the 
antenna  from  a  perfect >  reflecting  mirror.  We  need  to  account  for  the 
extra  dele*  and  the  change  in  direction  due  to  reflection.  Mart  again 
there  is  an  amplitude  factor  (the  sea*  one  that  is  ir-  (1))  which  Is  set 
to  unity: 

vricre  1  is  a  vector  direced  from  the  antenna  phase  center 
perpendicular  to  the  plane  of  the  reflector  and  having  a  length  equal  to 
the  closest  distance  of  the  reflector  plane  to  the  antenna.  tig  Is  the 
reflection  operator,  ar.  identity  transformation  that  projects  objects 
and  everything  else  into  'mirror  space': 

-si  -  ll"  ><  ^  \ 


The  ?-2»  factor  in  the  complex  exponential  accounts  for  the  excess 
propagation  delay.  We  need  tc  concern  ourselves  only  with 


I 
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suit! path  i roe  reflectors  placed  less  char,  a  code  chip  *way  in 
llgnt-tioe,  because  the  cross-correlatior.  properties  of  the  pseudorandom 
codes  automatically  reject  eultipath  froe  outside  a  chip  delay.  In  ocst 
cases,  in  fact,  we  need  to  worry  about  delays  of  only  a  snail  fraction 
of  a  chip  because  the  power  attentuatlon  on  the  smrts  of  the 
autocorrelation  function  of  the  code  is  relatively  high.  Also,  Equation 
(3)  neglects  dielectric  losses  of  nonideal  reflectors,  but  considers 
instead  onlv  perfect  reflection. 


The  total  received  signal  is  the  sub  of  Equations  (1)  and  (3):  A, ot 

•A  *  Aq.  For  the  general  case  of  a  readonly  oriented  and  shaped 
reflecting  surface  it  say  be  bo re  convenient  to  resort  tc  confutations! 
techniques  to  staplify  Ai0.  into  sagmtude  and  pnase.  Two  special 
cases  are  of  interest:  horizontal  and  vertical  reflect  me  planes. 
Generally  the  horizontal  caae  occurs  because  of  the  ground,  which  Is 
aostly  horizontal  and  so  nearly  sisilar  a:  the  two  baseline  cods  that 
the  effects  should  approxiaately  cancel  ir.  the  fornatioo  of  the 
observable;  thus  the  horizontal  reflector  case  is  fairly  benign  m  aost 
instances.  The  case  of  vertical  reflectors  is  often  aore  interesting 
and  less  benign  because  tbe  vertical  reflectors  are  not  necessarily 
present  at  both  antennas.  In  this  case,  for  a  reflector  m  the  plant  at 
distance  e  and  aziauth  t  the  reflector  orientation  vector  is 


e#v« 


O 
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an 6  substitution  into  Atct  gives  (after  ouch  sioplif icat log) : 

A*,,  *  ( '  * [  |  -  u*  f  «*?  (j  { ?Lm  -  4*)  -  *•« S <4-/*)})J 

•  Sin^'ti  *  4))-  k“-' 


The  urs  in  square  brackets  gives  the  deviation  (roe  the  expected 
antenna  pattern  Id.  Equation  (2)).  7hi»  ten  can  be  resolved  into 
polar  fora: 

|[*lj  *  »—**."*§  «•«*(-(/»■  “4)  '  b*n8  COb  (4  + 

(  b«C2(/a-4)  -  ^BST  ‘^i^ficoa/4"/^^)')  ^ 

4l-]  *  Un"'  f  t*.*1 bin  (*Q»-4)  ~  V^g  cgbCS^IjA 

l  -  W|c#s(20— 4)-  4^vn0cos(4-^))  J  • 

'  (5b) 


Tne  total  phase  angle  of  the  total  received  signal  is  the  phase  that  is 
ohser-ed  and  used  ir<  irterf erooetry .  This  total  phase  angle  is  found  by 
isolatinc  the  ergu&ents  of  the  cos? lex  exponential: 

^  %  *»»•  ^  W‘  j  Wj  \ 

**  *  v  j  “  605(14^-4)  - 

1:  :s  clear  that  the  pnase  error  is  *ust  the  expression  in  ($h,>.  This 
expression  has  beer,  plotted  in  Figure  j  fer  e*3C*  and  »*!.  and 

I'S  .  l **®rc .  with  <  varying  fro»  f*  tc  I8f*.  Tne  behavior  cf  the  phase 
error  ;s  oscillatory  for  a2*auths  well  a«ay  froo  opposition  (satellite, 

<»« 


ar.tenna,  and  cancer  of  tirror  copianar  in  vertical  plane).  One  can  see 
tn*  dependence  oi  cne  oscillation  :a:e  an t reaped  tc  aciouth  at  c  is 
variec:  cne  race  .*  higher  tor  oor*  cistart  re* lecicr*.  Ir.  the  region 
away  : ro=  opposition  this  error  free uency  say  be  high  enough  chat  the 
emulative  error  should  be  tero  s.  an  over  the  course  cf  an  observation 
session.  In  the  region  oi  opposition,  however,  the  phase  error  varies 
sore  slowly  w*tn  respect  it  ecioutr.  Ic  snou.d  be  nvte<  that  in  the 
course  of  as  hour  long  observation  a  CPS  sate. lice  typically  covers  onJv 
about  thirty  or  forty  degrees  ;n  annum.  It  is  possible  that  a 
satellite  couls  b«  observed  entirely  within  this  central  region  and  that 
the  sulcipath-re latec  error  could  contribute  a  non  cero  near,  phase 
offset,  in  the  next  section  w*  determine  when  this  observation  geonetry 
say  hold. 

Ceosetrv  of  hultipalfc 

Equation  '5b  gives  the  expression  for  phase  error  due  to  aultipair. 
given  chat  a  sate. life  located  al  »*.;  <**  be  observed  by  the  arter.r.a 

ir.  a  vertical  mirror  .ocated  at  annum  .  ar.c  distance  s.  «e  now  wish 
to  detenane  conditions  for  observability  sf  suit  .path  fur  this  specific 
reflecting  structure:  when  can  w*  see  the  satellite  ir  the  mirror' 

Scot  siople  ftoeetrit  constructions  will  show  when  this  situation 
arises.  First  recall  our  definitions  of  h  and  •:  is  the  atisutt  of 

the  vector  pointing  fr©«  the  antenna  phase  center  and  noma!  to  the 
plane  containing  the  reflector,  and  n  is  the  distance  fro*  the  antenna 
phase  center  to  the  piare  of  the  reflector  and  measured  along  the 
direction  cf  this  ncrsal  vector.  Our  observation  geometry  is 
iliustreted  in  Figure  *.  It  consists  of  a  vertical  mirror  of  height  f 
and  width  *-“>  fheight  is  oeasured  with  respect  ic  1  he  ground  plane 
.ocated  ir.  the  plane  described  V-  a  •  'ncos..  o*it*.  *  *  .  ^*e 

rectarguiar  nirmr  wiil  r»f.»ct  the  signal  frvt  *  setel.ite  at 
elevation  -  .  annum  ►  if  “tar.*- ;  ><t  ;  8f  r*an“  r  n  and  .<taf* 

'hers  ;  m  .  ■>«  viewing  rrctlen  c?  ideniita.  *.t  the  prti.et.  ff 

viewing  "f-  satriiite  fro*  the  isage  arter.r.a  ""be  hinc*  t**  n.  rrc  r 
through  a  w. id rw  located  where  the  nirrer  ;«  arc  rut  out  r.f  ar 

ooacue  dare. 


»  » 


A  staple  prifV.cal  t*chHi<ioe  tor  cetereir.ing  auiiipotr.  visibility 
suggest*  itselt.  tic  mg  a  polar  plot  cr  i>«uU:tr  I'cjrf.erjtr  in  the 
local  airy  tor  trie  rices  a:  writer  observations  are  »aae,  on*  can  ever. ay 
trie  windowing  tuner  ton  »  j  ;• ; » c  i&tt  pv.a*  cuorcinaies  tv  staving  a  curve, 
C,  ot  r  •  tan"  (nci>*(„  '  ®  tor  the  range  „  •  90*  <  »  '70*. 

Thu*.  it  tne  nirror  faces  north  ®lb<  ' *  ti*  curve  drawn  will  lie  or.  me 
north  halt  et  the  diagras.  A  satellite  a:  j*e*;iior;  tv.tsv'  t»r 
wtuch  <l#v*ttea  >ifv  it  lower  than  the  curve  C  ipesiuw.  outside 

:n*r  curve  C>  can  5*  *ee«  s#  a  vertic*.  reflector  of  height  h  and 
distance  s  fro*  the  antenna.  The  acisutha!  iirut*  tf  visibility  car,  he 
places  ha  drawing  ray*  at  s  *  *.„•;*?*  -  -  tar"  it  a  *n«  c  *  tw*4#0*}  * 
tan"*  i<*;  a, .  A  sat®. .it*  at  position  wrtcr  lie*  between  theta 

hound*  in  asisutn  and  outside  the  curve  C  can  h*  teen  m  the  plate 

*a.r? if  r  ;  C’-k  J  m  <k. 

\f»  e*,aPt>  ef  the  t  »;  JuSt  describee  4*  Sh<5*Vt»  in  figure 

*,  t ?a;*?t sf ••*  f^r  perinC#  r?  eiSlfellilv  uSeC 

t’f  l  fl  ■  i"  f  t  a  f  'lie  4  *  S  S  if.-.,  it  ,t  .  *  »»j  ;  r  ,  t.  fl  *  , , 

a  fli  *  ■' .  ‘ .  h  *  *  '  .  an*  •  *  -  ;  ) •  *  nfcirrer  fete*  northeast).  T1-® 
ahedui  ?dt(n  r  (it  •  t-  ttfcMtdb**  Ihee*  i* -*!■/»■.  t  *  JW*  - 

:  aft*'  ,  S  n  •  I  w  •  ta6M,?h  th  i  it  the  r  eg  ter  of  the  t*»  iP  which  a 

*<*:«.. it®  it  ntiiu*  i*  tsw  reflector.  In  thit  nashl*.  fvt  tt  the  only 
cat*  Hit*  that  frfca*®  error  due  tr  ovititHitr,  t®{  i*f  i  1  wt.s  fr«*« 

the  c*f  t .  -*i  plate. 

**vi ‘  1  r  a  :r  f^erer  1sh»tv* 

a  ?>.*«ri;oe  «nfv*t  ear® f ie*c t  *»«  f»*for»ei  t r  } ■'■'*-*«!  if  ***  th®  ®ff*ct#  cl 
i  -sn  1 1 1  pa?  *i  c>n  ;*>'  f*M«®  r*S,»e.  T^»«  *•«*■»■«'  *>•  tur  'vice, 

•■>«;•*  wkihPuS  a  fefl»c;iac  piste  and  '■w®  wit?  g  a  fwn.1  Mgh4  ;  f  r-<- 1  wide 

a.  icuirmi  piece  (>«;•«  f®et  ffofl  t‘«  art  «**•/* ,  t»e«5r.-®4  arid 

facing  aryfesiaetel®  s'*  w»«t  ef  n'Cth*  *  ;?**..  ?h*  c ta.hitai 

tec.hr -.sue  if  used  in  figure  *■  it  «» "n  th®  r®£irr  cf  *-.V»iaty  fry 

.rcw:  hur:e«  :n*  a;r*c-r  ‘we  ®cf*  »;  to®  *170'  a**.i®r ,r»* 

.ar®.  ; c r  L«  •  r®»;  r-  a  *,«®  >:i«  wa«  :cfl*is*r*^  i r  r*  » 

*a«*  r*f*ertcr  iia:.*eer.t. 


t  n*-v 


situation.  C*t  t .oe  satall  to*.  or.;y  $v<.  *nc  sva  are  visift^e  dunce  eoo* 
porttOB  of  ts«  osoervatioe.  p*r;oe.  arcaus*  JV*  ;s  at  a  iw«arT  «»ev anoc, 
we  *xp#et  it*  ux;wc  jma§e  e*vi*tier.  tt  S*  fc;£fc«*?  *t*ar.  t&at  o i  SV«, 

ftaaad  ee  is*  plot*  is  rtp*r*  >.  fc»tt  $v*«  par*  *.r.rou£f.  oppentic:.  «!r. 
(U  mirror,  wuu  to*  pna*#  arror  variac  with  aeivutr.  cooparee 

vritB  «i*a*n#r*.  W  S5CV.C  ;<*»*  tor  ifci*  *:#«.  v*rutnr  *!!*«  ;«<  t**« 
?eat«ftt  r«*ldw*.*.  for  fvrt*»#r  ’.w  iMte’.iai  ?*■•“  «m* 

'£<uaiies  }?’  tsr  :  tt*t  ft  «>•:*.;;:«  vj*.a;Jit>  .■  i*  piettoe  to? 

•(«  ttijtctftUf  of  SV >  aao  ,?  Fjp*r«  7R*  .arc®  e#viat Jv»  at  SV* 
ha*  a  »«a*  ie  :?owft  «*;»»*  »•  *?rr*>a;8at«i*  ;**  «r  *fc*»«t  ,m. 

iu»*i4a*:»  altar  a  iaaat  Sa’.tt  fit  to  p*a*a  •*a**r#a#«t*  are 

praaastatf  u  ftfwra  I  far  t*v  trace#*  (SV  t  . 

vttftaw:  is#  plat*  #8-  ;*  Cat  ♦  ,  a*a  astr.  taa  f.*l*  *»r.  ??  t*ct  %•»* 

i*  a  rroaowaeat  r*ru:u«  it  is*  f*  *;*»*:»  at  sv*  iiaa**  '.w#«r  t f*« 
Mttift  #tta  tA«  a*  rs.»  >*a«  t&  trowf*  «#*;#:. «•»  i* 

aa.  Si  so#  t»a  fit  iiftiaia*#  is#  #*t*#r*a  ®rr^>f  ft#  *c;**tiM  *•*«• 
5*a«;is#  »#tt  »r  ,.«#> ? s#s.tf  #»«f  #?  aaa;**  *  *i»*.  .:  *.f?i  ft*  »*?#<*♦* 
trtat  tfc#  #«..(  *#t  *?#«i**i*  tra«*  :*»*  t**#*r*t  .ca.  c«tr«. 

Ttvara  u  Mtan:  aa  apaarast;*  atfoef  tarraiaii#*  *«»#»<  it#  miew.t 
ot  Jt*  aa*  t a#  »r»cu:*<  error*.  ?a  tit*  at  I*#  r#*se#ai*  at  It*  tt» 
®**a  *ma.i#r  aa*  is#  e  #*&*>•  cti*«-  «Mft  !V  »r*<ici*i  #ft«t  •**  l#  t»ato*f 
t*  m*«.  n#  aia*  at  to#  r**’.«#*;«  f*r  •*«  »**  »*  «?  i**  »*»  «>**** 
«*  taal  af  Sv«  tft#*«  tv*  taaftet  ft*  ***^  it-  t*#  ftiat*.  Tt# 

u«i  ft#  a##  ;«  oi*#f  roofior  *ir#tt<*r**  *»#*#  a  c#ft«3aii#t  tat 
ft#  *e#r.  fra#  sr#  ills  t«  is#  ;**.►. 

ts«  *trui»#  *ft#«i  *t  :a#  m  t«»i«w;t  ;«  it*  *»*i  t*ari«i#f  ttcmeti t 
tor  at»r  a#*#;  at  7v«  s.ti  *f  f»«  ,#»a.*  miiitti**;* 

IS#  Sfa*«  »f*e*  fa.*.*  lie-eare*  TSaf*  #«*  otWt 

,*♦-*  *r*  s#t  a,r*»t  it#  tr#o 

?#n«-k3f  if  tr#  *r?;r  4xa*  ii  ;•*  fef  i‘i»  f#a*#t  ♦■f  *ot« 

ev»a-titati“*  arfwaests  s*  tr*s**t. 


##  -«v«  f *•  t '■*#■« ♦•ti i*  1  it* 


-.•at  *.t»  i*i* 


* 


!ui  tu  trrvr  w*  v  fiven  ta  t f rat  oj  tte  4eteu<i 

pattrfft  «.>e  t:«  ;.*iuor*  ot  i»-«  abc  tpeeal 

(4«t  u!  a  ;.«u  !<:.«((«(  «•»  ;a  toe*  etytt.  For  ;*;* 

6  f>«  »  i  - .  -  is  «  ,  4:  t4:  <J«f  Vr  «  «|*?  a  ii*  *  a£  A 

tuMtufiitt  it>  tdtMUBc  taee  im  trwuntei  ; » 

t.  u>«  if.  tu 

4»  *►*»  ..  Jvitti  *. .  .t.«»t.^tt  *«e  $;«♦  ef 

euitltAtS  !!»*  •  t«*tt}»*«  (tUtetet.  Th«? 

»Nu<*  A  tMAfCtUllljt  pfXllttX! 

V  4  *  >*•  »  . 

Tt>«  »>«’.  a*tu<  e«(  f  •*«*.<.  «4  tu  «<  S »  vjt*  t«t«ui 

Ji»w«e*»S  Wft4»f  hm  (4«  *11*  l  IM  *ai#MUM>  W*«<  W  fftaf 

*?  «at  **o*»r*. t i t 

ef  efr*  *«.*:n*»  ,»!«**  *n  »•(•<*  «»  *i*  }«,.«*  <t<* 

e«*ieu«  •«»».  ^tuf  **t*K>4  •*.*{*•  Nift  ;»♦**«  t(*t  »«■  U  te  <  t 

v.«  $&*»••  *?*••  **.<  ;t  ? ? »*  tu 

■*•?■>»«  j**  «■**?  *f*#f**  Muuii«*.  Ti-»*  uitx^  a«t*n4e* 

<■'*>  M#;»  #$  tu  t«;l9ttt*|  *«*?#<•♦ 

?eff#*s  ia»i*««t«t toe  ti  tu  mt*t«t  «u*t tu*  <  t a )  *«4  it) ) 

t»  ft  tu  »*«»•  *»iM  tu  t»  *  ttiMtui  utut*. 
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-  CIRCLED  NUMBERS  ARE  GPS  NAVSTAR  SATELLITE  DESIGNATORS  1*)  PM  Cfli£  m>**ber- 

-  SATELLITE  TRACKS  ARE  DARKENED  TO  INDICATE  SIMULTANEOUS 
VISIBILITY  ABOVE  20*  ELEVATION 


Figure  5.  Field  of  view  of  vertical  mirror  on  polar  plot  of  satellite 
positions  during  an  observation  session.  Analysis  is  for 
a  »  b  »  0.5,  a  ■  h  ■  1. 
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-  CIRCLED  NUMBERS  ARE  CPS  NAVSTAR  SATELLITE  DESIGNATORS  PRW calc  *H**bt«- 

-  SATELLITE  TRACKS  ARE  DARKENED  TO  INDICATE  SIMULTANEOUS 
VISIBILITY  ABOVE  20*  ELEVATION 


Figure  6.  Field  of  view  of  multipath  experiment  with  vertical  mirror 
placed  it  it  •  135*,  with  m  »  1.5,  h  ■  4,  a  -  b  »  2. 
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i'-as:  senates  phase  fit  for  oultipath  experiment 
w  tnoit  plate.  'b>  vith  plate.  Compare 
r  ve  »nd  SV®  *rith  predicted  error  *ho».*n  in 
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